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METHOD AND SYSTEM FOR SHORTENING CHANNEL IMPULSE RESPONSE 
USING TIME DOMAIN EQUALIZATION FILTER 

FIELD OF THE INVENTION 

The present invention relates generally to time domain equalizer filters (TEQ) 
5 used in digital communications systems, including discrete multi-tone (DMT), duplex 
systems, with equipment such as modems, customer premises equipment (CPEs), central 
office (CO) equipment, etc. The invention is used in physical layer applications and may 
be used in conjunction with other filters, such as echo cancellers. The invention is for use 
in many applications including with systems using asynmietric transmit and receive rates, 
10 e.g., ADSL, VDSL, etc., and in conjunction with a number of higher layer protocols, such 
as ATM, VoIP, VoDSL, etc. The present invention may further involve a minimum mean 
CI square error linearly constrained fast algorithm for adaptive training of time domain 

p equalizers. 

^ 15 BACKGROUND OF THE INVENTION 

si; 

K« The advent of the Internet and the widespread popularity of personal computers 

have created an unprecedented demand for high bandwidth networks. Generally, Internet 
5? applications, from simple email to real time video conferencing, from web surfing to 

U interactive movies, from interactive games to virtual TV stations, from online trading to 

20 online gambling, demand a higher bandwid± conununication network. A fundamental 
challenge for the conmiunication industry is to provide a reliable and affordable high 
bandwidth conmiunication link to all types of Internet users. Various competing wire- 
line, wireless, and optical broadband technologies are deployed to partially meet the ever- 
increasing demand for higher bandwidth. The fastest growing broadband technology is 
25 the Digital Subscriber Line (DSL) technology, which provides a high bandwidth always- 
on connection over standard twisted pair copper media of the conventional telephone 
network. Among other wire-line media, coaxial cables are capable of providing always- 
on connections, however, its presence is insignificant compared to millions and millions 
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of wired telephone customers who are connected by a twisted pair of copper wires. Other 
technologies, such as satellite, wireless, and optical, either provide limited coverage, 
limited bandwidth, or are too expensive for deployment to individual customers. As a 
result, DSL technology is uniquely positioned to provide the broadband link between 

individual customer premise and the central office, the so-called last-mile of the high- 
bandwidth conmiunication network, 

DSL is the fastest growing among emerging broadband technologies for very 
good reasons. First of all, DSL utilizes the existing copper wire network infrastructure. 
Secondly, compared to the voice modems, such as V.34 and V.90, used in most personal 
computers that provide up to 56 kbps dial-up connection, DSL provides a high bandwidth 
always-on connection with typical connection speeds from 384 kbps to 6 Mbps. 
Moreover, DSL is affordable with easy installation, simple turn-up, and high service 
reliability. The successful deployment of DSL is capable of providing digital broadband 
connection to anyone with an analog telephone line. 

DSL services have been standardized over time by regional organizations such as, 
American National Standard Listitute (ANSI), European Telecommunication Standard 
Institute (ETSI), and by world telecommunication organization International 
Telecommunication Union (ITU)- These DSL standards define data communication 
protocols to connect customer premise equipment (CPE) to the central office (CO) and to 
provide connections to various networks, such as DSL service providers, virtual private 
networks (VPN), or the Internet. Various forms of digital data (e.g., voice, video, and 
data) can be transported using DSL technology. For transport of voice, DSL equipment 
is connected to the public switched telephone network (PSTN). For transport of video 
and data, DSL equipment uses the Internet via an Memet service provider (ISP). Voice 
over DSL (VoDSL) is capable of providing computer-to-computer, computer-to- 
telephone, and telephone-to-telephone voice services using an integrated access device 
(IAD). Video over DSL includes transport of MPEG-1 or MPEG-2 files, video 
conferencing using Internet Protocol (IP) standard such as ITU H.323, WebCam, and 
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video mail. In addition, DSL supports simple data transport, e.g., bearer services, for 
virtual private network (VPN), leased data line such as Tl and El, Point-to-Point 
Protocol (PPP), Asynchronous transfer mode (ATM), and Internet Protocol (IP). 

Like other communication technologies, DSL has gone though a major evolution 
5 over the last decade and a collection of technologies, conmionly referred to as xDSL, are 
developed under the umbrella of DSL. One type of subscriber loop digital transmission 
technology involves an integrated services digital network (ISDN), which has replaced a 
significant portion of the analog phone lines in Europe and Japan. ISDN offers integrated 
voice and data services and connection speed up to 144 kbps. Due to the high cost of 
10 deployment, an alternative solution called integrated digital loop carrier (IDLC) was 
deployed in United States. However, resulting data rates were considered inadequate for 
fl^ individual customers. As a result, advanced DSL technologies were developed including 
S[ HDSL, SDSL, ADSL, HDSL2, SHDSL, and VDSL, all of which are capable of 

H; connection speed in excess of 1Mbps. These advanced DSL technologies were 

i=&!r. 

£ 15 developed to address different needs and application demands, while serving different 

J': market segments. For example, SHDSL is a symmetric service designed for long reach 

H= office applications with connection speed of 1.5Mbps, whereas, VDSL is designed to 

S provide a very high-speed asynmietric service for a short-range applications. 

DSL systems achieve high bit rates due to an effective appUcation of Digital 
20 Signal Processing (DSP) techniques which are enabled by advances in DSP algorithm 
design and VLSI design, which provides tools to physically implement these techniques 
and corresponding algorithms. Among the various techniques designed for 
communication applications. Discrete MultiTone (DMT) provides bit rates close to a 
maximum achievable level. In fact, due to this property and other properties, DMT has 
25 been elected as the line coding technique for the international ADSL standard. 

DMT technique deals with the linear distortion caused by channel spread where 
the channel is divided into smaller channels with no channel spread, thereby minimizing 
distortion. However, the length of allowable channel spread is limited by a symbol prefix 
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length used by the DMT system. In real life applications, however, the channel length is 
generally greater than the prefix length and further increases in length as the distance 
between the CO and the user increases. Therefore, a form of preprocessing received data 
is useful for shortening the effective channel length seen by the DMT receiver to a length 
5 less than or equal to the prefix length. The most common approach, known as Time 
domain Equalization (TEQ), is to apply a linear filtering to received data samples to 
shorten the effective channel spread. 

In a typical DMT based modem, the TEQ filter is located after the receive filter 
and before the DMT receiver block. Given this structure, one major concern involves 

10 designing an optimal and efficient algorithm to identify TEQ filter taps such that a 
requirement in the channel length is satisfied. Ultimate optimization goal involves 
maximizing data communication rate of the DMT system. However, designing TEQ 
filters to directly satisfy this goal is not practically implementable due to high 
computational complexity. Therefore, it is desirable to develop algorithms for obtaining 

15 TEQ coefficents which are lower in complexity and which provide a high bit rate. 



SUMMARY OF THE INVENTION 

Aspects of the present invention overcome the problems noted above, and realize 
additional advantages. One aspect of the invention provides a method and system to 

20 obtain Finite Impulse Response (FIR) filter coefficients for Time Domain Equalization 
(TEQ) for use in applications, such as Discrete Multi-tone (DMT) based applications, 
including asymmetric modes such as ADSL, VDSL, etc. The TEQ filter coefficients 
obtained by use of the invention shortens the overall effective discrete time channel 
impulse response length {hk » Zk) within a given target length (symbol/packet prefix 

25 length for DMT applications). Two key advantages, for example, of the invention are: 1) 
providing a TEQ filter with near optimal performance; and 2) facilitates adaptive 
implementation on a real system, an example of such an implementation is discussed 
elsewhere herein. Moreover, additional aspects of the invention include: a hypothetical 
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target impulse response (TIR) vector, representing the combined response of the 
transmission channel and TEQ filter, may be used as the search vector; one element, any 
element, of the TIR vector is constrained to be equal to a constant; the error is computed 
by using the difference between the TEQ output and the hypothetical TIR filter output 

5 and is minimized to lessen the residual energy of the overall channel outside the target 
length; TEQ filter coefficients may be computed by minimizing the mean square error 
(MSE) criterion; and a family of TEQ vectors may be parameterized by two quantities, 
the location of the constrained element of the TIR vector and the delay value used in the 
computation. The TEQ filter coefficients obtained as a result of the invention minimize 

10 energy ISI and ICI and yields a higher bit rate. 

In another manner, the present invention implements a minimum mean square 
error linearly constrained fast algorithm for adaptive training of time domain equalizers. 
One aspect of this invention relates to a minimum mean square error linearly constrained 
fast algorithm for adaptive training of a Time Domain Equalizer (MLC-TEQ) is 

15 provided. A fast adaptive algorithm of the present invention may be used to obtain Finite 
Impulse Response (FIR) filter coefficients for Time domain Equalizer (TEQ) used in 
Discrete Multitone (DMT) based applications, such as ADSL, for example. The TEQ 
coefficients obtained by the algorithm of the present invention shortens the overall 
effective discrete time channel impulse response length within a given target length (e.g., 

20 symbol prefix length for DMT application). Advantages of the proposed data aided 
adaptive algorithm may include providing the TEQ filter coefficients with near-optimal 
performance; having low computational requirements, having fast convergence, and 
exhibiting attractive stability properties. Other advantages may also be realized by the 
present invention and variations thereof. 

25 According to an embodiment of the present invention, a method for implementing 

a fast adaptive algorithm for obtaining finite impulse response filter coefficients for a 
time domain equalizer filter comprises the steps of adaptively computing at least one 
equalization delay parameter; and adaptively computing at least one time domain 
equalizer filter coefficient based on the equalization delay parameter. 
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Other aspects of the present invention may include an overall channel impulse 
response length shortened within a given target length; the step of computing an estimate 
cross-correlation function; the step of defining a variable as an argument maximizing an 
absolute value of the cross-correlation function wherein the variable represents a peak 
5 location for the absolute value of the cross-correlation function; the step of selecting an 
equalization delay as a function of the variable; the step of minimizing mean square error 
criterion; the step of implementing one or more of training sequences and received 
sequences; wherein the training sequences comprise consecutive samples of a received 
signal; wherein the time domain equalizer filter is a sample spaced finite impulse 
10 response filter; wherein the time domain equalizer filter is a fractionally spaced finite 
H; impulse response filter; and wherein the time domain equalizer filter minimizes one or 
p more of energy inter symbol interference and inter channel interference. 

m 

According to another embodiment of the present invention, a system for 
nl implementing a fast adaptive algorithm for obtaining finite impulse response filter 

J!: 15 coefficients for a time domain equalizer filter comprises a delay module for adaptively 
: computing at least one equalization delay parameter and an equalizer module for 

pd adaptively computing at least one time domain equalizer filter coefficient based on the 

Jyi equalization delay parameter. 

Other aspects of the present invention may include an overall channel impulse 
20 response length shortened within a given target length; a cross-correlation module for 
computing an estimate cross-correlation function; a defining module for defining a 
variable as an argument maximizing an absolute value of the cross-correlation function 
wherein the variable represents a peak location for the absolute value of the cross- 
correlation function; a selection module for selecting an equalization delay as a function 
25 of the variable; a minimizing module for minimizing mean square error criterion; an 
implementing module for implementing one or more of training sequences and received 
sequences; wherein the training sequences comprise consecutive samples of a received 
signal; wherein the time domain equalizer filter is a sample spaced finite impulse 
response filter; wherein the time domain equalizer filter is a fractionally spaced finite 
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impulse response filter; and wherein the time domain equalizer filter minimizes one or 
more of energy inter symbol interference and inter channel interference. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The present invention can be understood more completely by reading the 
following Detailed Description of the Invention, in conjunction with the accompanying 
drawings, in which: 

FIG. la is a block diagram illustrating a digital echo canceller, according to an 
embodiment of the present invention. 

FIG. lb is a system diagram illustrating a digital echo canceller implementation, 
according to an embodiment of the present invention. 

FIG. Ic is a system diagram illustrating a TEQ implementation, according to an 
embodiment of the present invention. 

FIG. Id is a system diagram illustrating modem transmission, according to an 
embodiment of the present invention. 

FIG. 2 is a block diagram illustrating an echo canceller, according to an 
embodiment of the present invention. 

FIG. 3 is a block diagram illustrating an example of a dual rate echo canceller, 
according to an embodiment of the present invention. 

FIG. 4 is a block diagram illustrating an example of a dual rate echo canceller, 
according to an embodiment of the present invention. 

FIG. 5 is a block diagram of a CPE side dual rate echo canceller for simultaneous 
support of multiple annexes, according to an embodiment of the present invention. 

FIG. 6 is an example of an impulse response of IF and AAF, according to an 
embodiment of the present invention. 
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FIG, 7 is a convolution of IF and AAF impulse response, according to an 
embodiment of the present invention. 

FIG. 8 is a block diagram illustrating an example of a weighted vector error echo 
canceller, according to an embodiment of the present invention. 

5 FIG. 9 is a flowchart illustrating an update process, according to an embodiment 

of the present invention. 

FIG. 10a is a block diagram of a TEQ system, according to an embodiment of the 
present invention illustrating a fractionally spaced TEQ. 

FIG. 10b is a block diagram of a TEQ system, according to an embodiment of the 
10 present invention illustrating a sample spaced TEQ. 

S FIG. 10c is a block diagram of a TEQ system, according to an embodiment of the 

O present invention illustrating a hypothetical model of a delay sample transmit signal for 

determining TEQ filter coefficients. 

S FIG. 11 presents a flow chart of the method in accordance with the current 

15 invention. 

^ FIG. 12 presents a flow chart describing a method of calculating the TEQ filter 

G vector of the present invention. 

FIG. 13 presents a flow chart of the steps of calculating TIR ({bk}) and TEQ 
({wk}) filter coefficients using MLC-TEQ. 

20 FIGs. 14a and 14b are schematic diagrams of hardware architectures in which the 

inventive aspects of the present invention may be incorporated. 

FIG. 15 is a block diagram of a physical media dependent layer of a ADSL CPE 
chip in which the inventive aspects of the present invention may be incorporated. 

FIG. 16 is a block diagram of a transmission convergence layer of a ADSL CPE 
25 chip in which the inventive aspects of the present invention may be incorporated. 
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FIG. 17 is a block diagram of an analog front end device in which the inventive 
aspects of the present invention may be incorporated, 

FIGs. 18a- 18c are block diagrams of applications in which the inventive aspects 
of the present invention may be incorporated. 

5 

LIST OF ACRONYMS 

A/D - Analog To Digital 
AAF - Anti-Aliasing Filter 
ADC - Analog to Digital Converter 
ADSL - Asymmetric Digital Subscriber Line 
AECF - Adaptive Echo Cancellation Filter 
AFE - Analog Front End 
AGC - Automatic Gain Control 
AOC - ADSL Overhead Control 
ANSI - American National Standard Institute 
ARM - Advanced RISC Machine 
ASIC - Application-Specific Integrated Circuit 
ATM - Asynchronous Transfer Mode 
BIU - Bus Interface Unit 
CIR - Channel Impulse Response 
CLECs - Competitive Local Exchange Carriers 
CLK - Clock 
CO - Central Office 
CPE - Customer Premise Equipment 
CRC - Cyclic Redundancy Check 
CRL - Clock Recovery Loop 
CSR - Control and Status Registers 
DAC - Digital to Analog Converter 
D/A - Digital To Analog 
dB - Decibels 

DMA - Direct Memory Access/Addressing 
DMT - Discrete Multi-Tone 
DP - DMT Processor 
DSB - Down Sampling Block 
DSL - Digital Subscriber Line 
DSLAM - DSL Access Multiplexor/Module 
DSP - Digital Signal Processing 
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• EEPROM - Electrically Erasable Progranmiable Read Only Memory 

• EC - Echo Canceller 

• ECF - Echo Cancellation Filter 

• EOC - Embedded Overhead Control 
5 • EPB - External Peripheral Bus 

• Eth PHY - Ethernet Physical Layer 

• ETSI - European Telecommunication Standard Institute 

• EWMP - Error Weighting Multi-input-multi-output Filter 

• FDD - Frequency Division Duplex 
10 • FEC - Forward Error Correction 

• FEQ - Frequency Domain Equalization 

• FIFO - First In First Out 

• FJR - Finite Impulse Response 

• GPIO - General Purpose Input/Output 
j:: 15 • HDLC-High-Level Data Link Control 

pi • HDSL - High-Speed Digital Subscriber Line 

• HIU - Host Interface Unit 

O • IAD - Integrated Access Device 

m • ICE - Information and Content Exchange 

20 • ICE -In-Circuit Emulator 

• ICI - Inter Channel or Carrier Interference 

• IDLC - Integrated Digital Loop Carrier 
m •IF - Interpolation Filter 

H • VF - Interface 

IB 25 • IFB - Interpolation Filter Bank 

Q • IFFT - Inverse Fast Fourier Transform 

^ • ILECs - Incumbent Local Exchange Carriers 

• IP - Internet Protocol 

• IPF - Interpolation Filter 

30 • ISDN - Integrated Services Digital Network 

• ISI - Inter Symbol Interference 

• ISOS - Integrated Software on Silicon 

• ISP - Internet Service Provider 

• rrU - International Telecommunications Union 
35 • DCCs - Interexchange Carriers 

• JEDEC - Joint Electronic Device Engineering Counsel 

• JTAG - Joint Test Action Group 

• LD - Line Driver 

• LMS - Least Mean Square 
40 • LPF - Low Pass Filter 

• LTI - Linear Time-Invariant 
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• MAC - Media Access Control 

• MCM - Multi Chip Module 

• Mil - Media Independent Interface 

• MIPS - Million Instructions Per Second 
5 • MIMO - Multi-Input Multi-Output 

• MMSE - Minimum Mean Squared Error 

• MSE - Mean Squared Error 

• NP - Network Processor 

• NTR - Network Timing Reference 

10 • 0AM - Operation, Administration and Maintenance 

• OSI - Open Source Initiative 

• PAC - Programmable Attenuation Control 

• PCI - Peripheral Component Interconnect 
. . • PGA - Progranmiable Gain Amplifier 

pj 15 • PHY -Physical Layer Device 

3 • PMS - Physical Media Specific 

fU • PP - Protocol Processor 

CI • PPP - Point-to-Point Protocol 

n.- • PROM - Programmable Read Only Memory 

20 • PSD - Power Spectral Density 

f" • PSTN - Public Switched Telephone Network 

• RCDR-EC - Reduced Complexity Dual Rate Echo Canceller 

• RMB - Rate Matching Block 

M: • RMB-RT - Rate Matching Block from Receive to Transnait 

W 25 • RMB-TR - Rate Matching Block from Transmit to Receive 

p • ROM - Read Only Memory 

• RS - Reed Solomon 

• RSIF - Reference Signal Interpolation Filter 

• Rx - Receive 

30 • SD - Sigma-Delta 

• SDRAM - Synchronous Dynamic Random Access Memory 

• SNMP - Simple Network Management Protocol 

• SNR - Signal To Noise Ratio 

• SRAM - Static Random Access Memory 
35 • STM - Synchronous Transfer Mode 

• S VD - Singular Value Decomposition 

• TC - Transmission Convergence 

• TEQ - Time Domain Equalizers 

• TIR - Target Impulse Response 
40 • Tx - Transmit 

• UART - Universal Asynchronous Receiver Transmitter 
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• USB - Up Sampling Block 

• USB - Universal Serial Bus 

• VCXO - Voltage Controlled Crystal Oscillator 

• VDSL - Very High Bit Rate Digital Subscriber Line 
5 • VoDSL - Voice over Digital Subscriber Line 

• VoIP - Voice over Internet Protocol 

• VPN - Virtual Private Network 

• VU - Vectorization Unit 

• WEVE-EC - Weighted Vector Error Echo Canceller 

10 

DETAILED DESCRIPTION OF THE INVENTION 

The following description is intended to convey a thorough understanding of the 
invention by providing a number of specific embodiments and details involving echo 

15 cancellers. It is understood, however, that the invention is not limited to these specific 
embodiments and details, which are exemplary only. It is further understood that one 
possessing ordinary skill in the art, in light of known systems and methods, would 
appreciate the use of the invention for its intended purposes and benefits in any number 
of altemative embodiments, depending upon specific design and other needs. 

20 According to an embodiment of the present invention, a dual rate echo canceller is 

provided for applications with asymmetric transmit and receive rates. In particular, the 
present invention provides Reduced Complexity Dual Rate Echo Canceller (RCDR-EC) 
architecture along with low complexity Least Mean Square (LMS) update rules. 
According to one aspect, the present invention is directed to an echo canceller that 

25 provides rate matching functionality for applications (e.g., ADSL, VDSL, etc.) with 
asymmetric data rates. 

The RCDR-EC architecture and LMS update rules of the present invention 
provide various features and advantages. For example, RCDR-EC may operate at the 
lower of the two rates, receive and transmit, requiring less computation per data sample. 

30 RCDR-EC implements a significantly smaller length echo cancellation filter (ECF) for 
achieving the same (or similar) level of echo suppression of conventional 
implementations. This reduces the hardware requurement for implementation of RCDR- 
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EC. Due, in part, to the reduction in ECF length, RCDR-EC may involve significantly 
less computation for filtering operation in steady state. The reduced rate and smaller 
echo canceller length enables RCDR-EC training to be simplified where less computation 
is required for adaptive training of the RCDR-EC. In addition, inherent out-of-band 
5 noise rejection capability of RCDR-EC enables adaptive training to require less time to 
converge. Moreover, out-of-band noise rejection enables the adaptive training algorithm 
to yield an improved RCDR-ECF resulting in enhanced echo suppression. 

The present invention provides an efficient rate matching echo canceller for 
applications with asymmetric transmit and receive rates, such as ADSL and VDSL, for 
10 example. Another aspect of the present invention provides echo canceller structures for 
cases involving higher and lower receive rates as compared to the transmit rate. Another 
aspect of the present invention is directed to developing a low complexity LMS update 
p for adaptive training of the echo cancellation filters. In addition, the present invention 

Ij reduces hardware requirements for implementation of the echo canceller and reduces 

=1== 15 computational requirements for steady state operation of the echo canceller. As a result, 
y, overall performance of the echo canceller is improved. 

A basic block diagram of the digital echo canceller structure 100, which may be 
m employed in central office (CO), CPE or other equipment, is shown in FIG. la. A 

rr transmit filter 110 receives an input transmit signal from a transmit path and generates a 

20 filtered transmit signal The output of transmit filter 110 may be received by analog front 
end (AFE) 112, which is comprised of digital to analog converter and various analog 
components including line drivers, for example. Echo cancellation filter (ECF) 122 
receives an input from the output of the transmit filter and generates a copy of an echo 
signal by a linear filtering operation. The output generated by AFE 112 may be received 
25 at an input of Analog Hybrid 114, which is comprised of various passive analog 
components, for example, and may be coupled to a twisted wire pair or other 
transmission line. A receive signal may be received by analog hybrid 1 14, which may be 
coupled to AFE 116 at a receive side. For example, analog hybrid 114 may include 
common analog components, such as resistors and capacitors. The output generated by 
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AFE 116 is received at an input of receive filter 118. The receive filter perfomis linear 
filtering operation and attempts to suppress out of band noise. The copy of the echo 
signal generated by ECF 122 is subtracted from the receive signal generated by receive 
filter 118, at summer 120. The resulting residual signal may be used to train echo 
5 cancellation filter 122, such as by feeding a sampling signal 121 into echo cancellation 
filter 122. 

As shown in FIG. la, digital echo cancellation may be achieved by obtaining a 
reliable copy of the echo signal, which may be generated from the transmit signal using a 
filter, such as a Finite Impulse Response (FIR) ECF, as illustrated by 122. The copy of 
10 the echo signal is then subtracted from the received signal. The FIR ECF 122 operates as 
a connecting branch between the transmit and receive paths. For example, the echo 
canceller filter receives an input from the transmit path and generates an output for the 
receive path. Hence, the ECF 122 accommodates the rate transition from the input 
pj transmit rate to the output receive rate. For applications, such as symmetric DSL (e.g., 

^1, 15 HDSL2, G.shdsl, etc.), where the transmit and receive data rates are identical there is 
^ generally little or no need for rate transition. As a result, implementation of ECF is 

■^1 straightforward in these cases. However, for applications such as ADSL and VDSL and 

^ for other applications that support asymmetric transmit and receive data rates, rate 

□ transition may be required at the ECF, As addressed by the present invention, there are 

20 different approaches to achieve the rate transition within the ECF. In accordance with an 
embodiment of the present invention, ECF 122 may incorporate the inventive features as 
discussed in detail below. 

FIG. lb is a system diagram illustrating a digital echo canceller implementation, 
according to an embodiment of the present invention. In particular, system 130 illustrates 
25 echo cancellation functionality in a DSL modem system. CO side DSL modem is 
represented by 132, which may include a DSL modem 134 transmitting to and receiving 
from Hybrid and Analog Front End 138. Echo cancellation filter 136 compensates for a 
portion of the transmit signal that is received as an echo signal as shown by arrow 137. 
CO side DSL modem 132 may be connected to CPE DSL modem 142 by a twisted 
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copper pair 140. CPE DSL modem 142 may include a DSL modem 144 transmitting to 
and receiving from Hybrid and Analog Front End 148. Echo cancellation filter 146 
compensates for a portion of the transmit signal that is received as an echo signal as 
shown by arrow 147. The inventive aspects of the echo cancellation filter as described in 
detail below may be incorporated into echo cancellation filter represented by 136 and 
146. 

FIG. Ic is a system diagram illustrating a TEQ implementation, according to an 
embodiment of the present invention. In particular, system 150 represents TEQ in a 
DMT based DSL modem. CO Side DSL Modem 152 may include DMT transmitter 154, 
Transmit Filter 156, Analog Front End 158, Receive filter 160, TEQ 162 and DMT 
receiver 164. CO Side DSL Modem 152 is coupled to CPE DSL modem 168 via twisted 
copper pair 166. CPE DSL modem 168 may include DMT transmitter 170, Transmit 
Filter 172, Analog Front End 174, Receive filter 176, TEQ 178 and DMT receiver 179. 
The inventive aspects of the TEQ features as described in detail below may be 
incorporated into TEQ represented by 162 and 178. 

FIG. Id is a system diagram illustrating modem transmission, according to an 
embodiment of the present invention. At the CO end, various sources of data may be 
transmitted via a DSL modem, as represented by 182, which may further include DSL 
Access Module, DSLAM, for example. Intemet 184 may provide data, such as video and 
other forms of data, via Intemet Service Provider 183. ATM Network 186 may provide 
data, such as video and other forms of data, via ATM Switch 185. Public Switch 
Network (PSTN) 188 may provide voice data to Central Office Switch 187, DSL modem 
182 may receive the various forms of data and communicate via a twisted copper pair 
181 at a CPE device. At the CPE end, data from DSL modem 182 may be received by 
DSL modem CPE 192. Voice data may be received by splitter 190 for transnGiission to a 
telephone 191 or other voice device. DSL modem CPE 192 may transmit data to a 
computer 193 or other device. Video data may be received by set top box 194 for 
transmission to television (TV) 195 or other video device. The inventive aspects of the 
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echo cancellation features as well as the TEQ features may be incorporated in this 
system, in accordance with the various embodiments of the present invention as 
discussed in detail below. 

Generally, the receive data rate is a rational multiple of the transmit data rate, as 
illustrated by 

where a = P/Q for co-prime integers P and Q (e.g., P and Q do not have any 
common factors). Rtx and Rrx represent the sampling rates at the output of the transmit 
and the receive filters, respectively. In one example, a > 1 may represent cases, such as 
ADSL CPE, where the receive rate is higher than the transmit rate. On the other hand, a 
< 1 may represent cases, such as ADSL CO, where the receive rate is lower than the 
transmit rate. Depending on the value of a, different Reduced Complexity Dual Rate 
Echo Canceller (RCDR-EC) structures may be implemented. 

FIG. 2 is a block diagram of an echo canceller, according to an embodiment of the 
present invention. Sampling may occur by rate matching as represented by 220. In one 
example, FIG. 2 represents an echo canceller for applications where the receive rate is 
higher than the transmit rate. A transmit filter 210 may receive a transmit signal and 
generate a filtered transmit signal. The output of the transmit filter 210 may be coupled 
to an analog front end or other device. A copy of an echo signal associated with the 
transmit signal may be up-sampled by Up Sampling block 212 to match the sample rate 
to the receive rate. The up-sampled signal may then be received by Echo Cancellation 
Filter 214 to produce an output at the receive rate. The resulting output may be 
subtracted from a received signal filtered by Receive filter 216. Conventional echo 
cancellers generally operate the ECF at the higher receive rate. As a consequence, the 
ECF, which may be a Finite Impulse Response (FIR), may require a larger number of 
taps to span a fixed length of time. As a result, extra hardware or software complexities 
will result. In accordance with the present invention. Echo Cancellation Filter 214 may 
include the inventive aspects of the invention, as discussed in detail below. 
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FIG. 3 illustrates a block diagram of a Reduced Complexity Dual Rate Echo 
Canceller (RCDR-EC), according to an embodiment of the present invention. In 
particular, FIG. 3 illustrates a case where the receive rate is higher than the transmit rate. 
As shown in FIG. 3, a transmit filter 310 may receive a transmit signal from a transmit 

5 path and generate a filtered transmit signal x{n). Transmit filter 310 may operate at a 
transmit rate (Rtx)^ The output x{n) of the transmit filter 210 may be received by an 
analog front end or other device. Echo Cancellation Filter (ECF) 312 may include an 
adaptive FIR filter that receives an input from the transmit path and generates a copy of 
an echo associated with the transmit signal, as shown by xl(ny ECF 312 may operate at 

10 transmit rate (Rtx), Since both the input and output of ECF 312 are sampled at the lower 
transmit rate (Rtx), the ECF 312 operates at the lower transmit rate (RtxX 

Rate Matching Block from Transmit to Receive (RMB-TR) 320 up-samples the 
output xl{n) of ECF 312 by a factor a with appropriate filtering wherein a = P/Q. The 
RMB-TR 320 may fiirther include sub blocks, such as an Up Sampling Block (USB-P) 

15 322 with up sampling factor P, an Interpolation Filter (IFF) 324, and a Down Sampling 
Block 326 (DSB-Q) with down sampling factor Q. Other sub-blocks may be 
implemented. 

An output xl{n) of the ECF 312 may be received by the USB-P 322. An up 

sampling function at the USB-P 322 may be achieved by a zero filling operation, e.g., the 
20 USB-P inserts P-l zeros between consecutive input signal samples. Other operations for 

achieving up sampling may be implemented. Thus, USB-P 322 may generate samples P 

times faster than the input sample rate Rtx* 

An output x2(Z) of the USB-P 322 may be received by the IPF 324. The IPF 324 

may include a fixed FIR filter that interpolates the up-sampled signal. Essentially, IPF 
25 324 may serve as a low pass filter for filtering out most or all P-l replicas of the transmit 

signal generated from the up sampling function. The IPF 324 may operate at P times the 

transmit rate Ru or other variation thereof. 
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An output x3{l) of the IPF 324 may be received by the DSB-Q 326. The DSB-Q 
326 may down sample the output of the IPF 324 by a factor Q by dropping Q-l samples 
from the input data stream and generating a signal at the receive rate Rrx^ 

An incoming receive signal from an Analog Front End may be received by 
5 Receive Filter 314, which operates at a receive rate Rrx- Receive Filter 314 may then 
generate a filtered receive signal y(m). The output x4(m) of the RMB-TR 320 may then 
be subtracted from the filtered receive signal y(m) from Receive Filter 314 at sunmier 316 
to generate a residue echo signal erim). 

Rate Matching Block from Receive to Transmit (RMB-RT) 330 may down 
10 sample the residue echo signal erim) by a factor a (or other factor) with appropriate 
filtering. Similar to RMB-TR 320, the RMB-RT 330 may include sub blocks, such as Up 
Sampling Block (USB-Q) 332 with up sampling factor Q, an Anti-Aliasing Filter (AAF) 
334, and a Down Sampling Block (DSB-P) 336 with down sampling factor P. 

The residue echo signal er(m) may be received by the USB-Q 332. The up 
15 sampling at the USB-Q 332 may be achieved by a zero filling operation, e.g., the USB-Q 
332 inserts Q-l zeros between consecutive input signal samples. Other operations for 
achieving up sampling may be implemented. Thus, USB-Q 332 may generate samples Q 
times faster than the input sample rate Rrx* 

An output el(l) of the USB-Q 332 may be received by the AAF 334. The AAF 
20 334 may include a fixed low pass FIR filter that receives the up-sampled output of the 
summer 316 and filters out most or all the high frequency signals above the transmit band 
before feeding it to a down-sampling block. The AAF 334 may operate at Q times the 
receive rate Rrx^ 

An output 62(1) of the AAF 334 may be received by the DSB-P 336. The DSB-P 
25 336 may down sample the output of the AAF 334 by a factor P (or other factor) by 
dropping P-1 samples from the input data stream and generating an error signal e{n) at 
the transmit rate Rtx^ 
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In one particular embodiment, for RCDR-EC of the present invention, the up 
sampling is accomplished after the ECF 312, as opposed to the conventional approach 
where up sampling is done before EEC 312. 

The transmit signal x{n) may first be filtered by a N-tap FIR ECF where 

represents N adaptable ECF filter coefficients. The output of the ECF may be written as 

N 

xl{n) = ^ X (fi - A: + l)wfc = x\w, 

where Xn is the data vector 

Xn [xin\ xin - 1), x{n - 2), . . . jc(n - A/' + 1)]^, 
p« JO w is the coefficient vector, 

|5 and T denotes the transpose operation. 

According to another embodiment of the present invention, the RMB-TR 320 may 
up-sample xi(n) by a factor a. For example, for every P input samples the RMB-TR 320 
15 may generate Q output samples where a = P/Q. As discussed above, the up sampling at 
the USB-P 322 may be achieved by a zero filling operation. Other operations for 
achieving up sampling may be implemented. The output x2{l) of the USB-P 322 may be 
expressed as 

ifl = Pn 



ly 

Ms 



20 



\x\{n\ 



x2{t) = 

otherwise 



In other words, x2([) may include the samples of xl{n) with P - 1 zeros inserted in 
between, e.g., samples of x2([) may be given as follows: 

•• •0,jcl(n - 1),0,0,- • -,0, A;l(n),0,- • • 

V ' 

P-1 zeros 

Different time indices may be used to represent different sampling rates, e.g., n, /, 
25 and m may be used to represent rates Rtx, PRtx (or equivalently, QRrx\ and Rrx, 
respectively, for example. 

-19- 



PATENT 

Attorney Docket No. 56162.000322 



The output of a A^/tap FIR IPF 324 may be given by 

x3il)= J^x2(l-k+m, 

it=l 

where 

represent the Nf IPF coefficients. 

The output x3(l) of the IPF 324 may be down sampled using DSB-Q 326 where 
the output may be represented by 

x4(m) = x3(0, where / = Qm, 

In other words, x4(m) may include decimated samples of x3(n), where the samples 
of x4{n) may be written as 

.^.,x3((2(/- l)),x3(2/),x3(G(/+ 1)), ... 
The output of RMB-TR 320 may be subtracted from the incoming receive signal 
y(m) from 314 to generate the residue echo signal e^m), which may be defined as 

er{m) = y(m) - x4(m). 
The residue echo signal er(m) may be up-sampled by USB-Q 332 by zero filling 
(or other operation) where the output may be given by 

0 otherwise. 



The up-sampled signal el(/) may be filtered by the AAF 334, which may be a FIR 
filter, to generate the filtered version e2(Z), which may be expressed as 

e2(/)= £el(/-fe + l)g„ 

k=l 

where 

{gug2. g3, ... gNg} 

are a Ng AAF coefficients. 

The output e2(l) of the AAF 334 may be down sampled using DSB-P 336 to 
generate the error signal e(n) which may be defined as 
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e{n) = el{l\ where / = Pn. 

IPF 324 and/or AAF 334 may include low pass filters that filter out the high 
frequency signals above the transmit band. According to an embodiment of the present 
invention, the same (similar or related) coefficients may be used for both IPF and AAF, 
^•g-» fk = gk for all k. Also, the samples of the error signal may be generated at the 
transmit rate Rtx. The error signal e{n) may then be used to adaptively train the 
coefficients of ECF 312, as discussed below. 

To describe the LMS update rules for the RCDR-EC structure of the present 
invention, the error signal e{n) may be expressed in terms of the transmit signal x{n) and 
the receive signal y{m) as well as IPF and AAF coefficients. For example, the error 
signal e{n) may be written as 

^(«) = yr(n)-h^X„w, 

where yr(n) is the received signal component after USB-Q 332, AAF 334, and DSB-P 
336, w is the ECF coefficient vector, 

x{n-N + \) 
x{n-N) 

x(n-L-N + 2) 

is the L X N input data matrix, and 

h [huh2,h, ... htf, 
is the filter coefficients vector that represents the combined effect of IPF 324 and AAF 
334 along with the effects of USB-P 322, DSB-Q 326, USB-Q 332, and DSB-P 336. 

The DSB-Q 326 and USB-Q 332 may have a sub-sampling effect on the IPF 324 
and AAF 334 coefficients, e.g., the equivalent filters may include one of the possible Q 



x{n) x(n-l) 
x(n - 1) x{n - 2) • • • 

x{n-L + l) x(n-L) 
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phases of the original filters. For example, /^^ may denote one of possible Q phases of 
the IPF filter 324, where 

fk^^fQ(k-iMu forA:=l,2,...iV/, 

In addition, q\ % \, 2, 3, , Q denotes the selected phase where is the smallest 
integer not less than Nf/Q. Similarly, for the AAF coefficients, 

gk^ = gQ(k-i)+qU forA:=: 1,2, ...A/'/ 

where ^2 X 1, 2, 3, ... , g denotes the selected phase and A^^^ is the smallest integer not 
less than N/Q. 

The (A^ + Ng^ - 1) coefficients of the combined filter response of the single 
phase IPF 324 and AAF 334 may be represented by bi. For example, bi may be obtained 
by convolving/^^ and gp as follows: 

bi = X/*^ ^ = 1' 2, N^f-^lsfig- 1. 

k 

From the combined filter coefficients Z?/, the elements of h may be easily obtained. 
The effect of USB-P 322 and DSB-P 336 may include the selection of a phase out of 
possible P phases. In other words, 

hi = bpQ.x)^p, for /= 1,2, ... ,L, 

where p x 1, 2, 3, ... , P denotes the selected phase and L is the smallest integer not less 
than {Np- + Ng^ - 1)1 P. Though the exact value of /?, qu and ^2 may not significantly 
impact the performance, the a priori knowledge of these integers may be used to 
calculate h. In the equations, the value in parenthesis generally refers to the time and 
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variables in bold generally refer to arrays (e.g., vec(n) refers to values of the array "vec" 
at time n). 

Once the expression for the error signal e(n) is obtained, the LMS update rule 
may be given as follows: 

Wn+i =w„-/^e(n)X/h, 

where /u is the LMS step size. Time subscript for Wn may represent the LMS iteration. 

At a first glance, it may appear that the proposed LMS update requires (L + l)Ar 
multiplications per iteration where the factor L is due to the vector-matrix product X„^h. 
However, due to the shift structure of X„, e.g., X„ may be a Hankel (or other) matrix, per 
iteration, N multiplications may be used instead. Hence, the over all computational 
complexity of the LMS for RCDR-EC may be 2N multiplications per iteration. The 
multiplication required may be twice that of usual LMS which in comparison requires N 
multiplications. As a result, there is a savings of filter coefficients by, at least, a factor a. 
Thus, the computational complexity of the proposed LMS for RCDR-EC of the present 
invention is lower for a > 2. 

According to another embodiment of the present invention, a simplified LMS may 
be implemented to reduce computational complexities by using an approximation on h. 
In one example, IPF 324 and/or AAF 334 may include low pass filters. Hence, the 
impulse responses may be similar to a sine function. In other words, the impulse 
response has a dominant peak and diminishes in magnitude away from the peak. 
Moreover, the convolution of two such filters may also exhibit a similar behavior. 
Variations may be made to h. For example, the largest element of h may be replaced 
with its sign. In another example, some or the rest of the elements may be set to zero (or 
other value). Replacing the largest magnitude element by its sign may result in a scaling 
of h where such scaling may be absorbed in In addition, no approximation may be 
involved in this step. After such scaling, the remaining elements may have a magnitude 
less than one, which may then be approximated to zero. With this approximation, the 
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coefficient vector h becomes a vector of zeros except one non-zero element (e.g., ± 1). 
Without loss of generality, the largest magnitude element of h may be positive and the 
matrix-vector product may reduce to a delayed version of the input vector, e.g., 

x(n - d) 
x{n -d-l) 

x{n-d-N + \) 
/\ 

where h is the approximate h and d is the index of the largest element of h. 
Using this approximation, the LMS update may simplify to, 

JT- 10 w„+i =w„-/fe(n)x„.<i, 

r;^ where Xn is the data vector defined before as 

U Xn = [x(n\ x(n - lXx(n -2), ... x(n- N + 

5 Hence, the simplified LMS update for RCDR-EC may be similar to the normal 

LMS update except for the delay. In other words, the computational complexity of the 
simplified LMS for RCDR-EC is similar to that of usual LMS. 

Since both IPF and AAF are part of the system design, the approximation does 
20 not restrict the use of the simplified LMS. The IPF and AAF may be defined such that 
the resulting combined filter has a single dominant coefficient. Factors that may degrade 
the convergence and performance of LMS may include multiple dominant peaks in the 
combined filter response and inexact setting of the delay d, for example. 

Since up sampling in the USB-P may be performed by zero filling, the sub blocks 
25 of RMB-TR 330 may be combined for efficient poly-phase implementation. The poly- 
phase implementation of RMB-TR may require RrxN/P (or equivalently, RtxN/Q) 



5 Xn^h = 
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multiplications per second. Efficient poly-phase structure may also be used for RMB-RT 
330 and the number of multiplications that may be required per second is RtxN/Q. 

Various comparisons may be made between a conventional echo canceller and the 
RCDR-EC of the present invention for the case a > 1, for example. 
5 As for structural differences, a conventional echo canceller is generally composed 

of an up sampling block followed by an adaptive FIR ECF that operates at the higher 
receive rate Rj^. In contrast, in the RCDR-EC approach of the present invention the 
adaptive FIR ECF may operate at the lower transmit rate Rtx> A rate matching block 
RMB-TR 330 may be used to generate echo samples at the receive rate and another rate 
10 matching RMB-RT 320 may be used to calculated the error signal at the transmit rate. 

As for filter lengths, to capture the same (or similar) fixed amount of time span 
Q (e.g., echo tail), the RCDR-EC of the present invention may use filter taps as compared 

5 5 

pj to the oN taps required by the conventional echo canceller approach. However, in 

addition to the ECF, RCDR-EC of the present invention may implement additional filters, 
H 15 which may include DPF and AAF. In general, the ECF is longer than the combined length 
J'' of IFF and AAF. Hence, there is a significant reduction in the number of filter taps. For 

H; example, for a typical value of N = 256 and a = 8, the conventional filter requires 

approximately 2048 taps. According to an example of the present invention, even with 
^ two 64 taps DPF and AAF filters, the number of taps for RCDR-EC is approximately 384, 

H 20 significantly less than that required by conventional systems. 

As for LMS training computations, the computational complexity may be 
approximately equal (or similar) for the conventional LMS and simplified LMS for 
RCDR-EC of the present invention. However, since the conventional ECF is a times 
longer than that of RCDR-EC of the present invention, the number of multiplications 
25 required for LMS update may be a times higher for conventional echo canceller. Due to 
the reduction in the filter taps, LMS for RCDR-EC of the present invention has lower 
complexity then the conventional approach for a > 2, according to one example of the 
present invention. 
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As for steady state computational complexity, in steady state using efficient poly- 
phase implementation the total number of multiplications per second that may be required 
by RCDR-EC is (A^ + N/Q + N/Q)Rtx. In contrast, even with efficient poly-phase 
implementation, the number of multiplications per second required by the conventional 
5 echo canceller is at least oNRrx. Li a typical application, the total filter length of IPF and 
AAF may be considered small compared to the length of the ECF filter, e.g., NfhNg< N, 
Hence, RCDR-EC of the present invention may reduce the steady state computation by at 
least a factor of a^. 

As for improved noise floor level and overall performance, in RCDR-EC the AAF 
10 filters out signals above a transmit band before generating the error signal e(n). Thus, the 
error signal has a reduced contribution from the wide band receive signal which acts as 

Q noise for echo cancellers, hi a conventional echo canceller, the receive signal, even the 

D 

i^i part outside of the transmit band, contributes to the error signal. Hence, RCDR-EC of the 

y present invention has a lower noise floor that results in reduced convergence time and 

M. 15 improved overall performance of the LMS training algorithm. Other comparisons as well 
y as advantages may be recognized by the present invention and variations thereof. 

FIG. 4 is a block diagram of a Reduced Complexity Dual Rate Echo Canceller 
Lt (RCDR-EC), according to an embodiment of the present invention. The block diagram 

of FIG. 4 may support applications, such as ADSL CO, where the receive rate is lower 
20 than the transmit rate, e.g., a < 1. 

Transmit Filter 410 may receive a transmit signal from a transmit path and 
generate a filtered transmit signal x(n). The filtered transmit signal x(n) may be an input 
to an Analog Front End or other device. 

Rate Matching Block (RMB) 420 may up-sample the input signal by a factor a (or 
25 other factor) with appropriate filtering. The RMB 420 may include sub blocks, such as 
an Up Sampling Block with up samphng factor P (USB-P) 422, an Anti- Aliasing Filter 
(AAF) 424, and a Down Sampling Block with down sampling factor Q (DSB-Q) 426. 

The filtered transmit signal x(n) may be received by the USB-P 422. The up 
sampling at the USB-P 422 may be achieved by a zero fiUing operation, e.g., the USB-P 
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422 may insert P-l zeros between consecutive input signal samples. Other operations for 
achieving up sampling may be implemented. Thus, USB-P 422 generates samples P 
times faster than the input sample transmit rate Rtx. 

An output xl(n) from the USB-P 422 may serve as an input to the AAF 424. The 
AAF 424 may include a fixed low pass FIR filter that receives the up-sampled signal 
xl(n) and filters out most or all the high frequency signals above the receive band before 
feeding it to a down-sampling block. AAF 424 may operate at P times the transmit rate 

An output x2(n) from the AAF 424 may serve as an input to the DSB-Q 426. The 
DSB-Q 426 may down sample the output of the AAF 424 by a factor Q by dropping Q - 
1 samples from the input data stream and generate the output signal z(m) at the receive 
rate i?^. 

Echo Cancellation Filter (ECF) 428 may include an adaptive FIR filter that 
receives the down sampled signal z(m) and generates a copy of the echo. Since both the 
input and output of the ECF 428 are sampled at the lower receive rate (/?^), the ECF 428 
operates at the receive rate Rrx. According to an embodiment of the present invention, 
RCDR-EC provides down sampling before the ECF, wherein ECF operates at the lower 
receive rate 

The transmit signal x(n) may be up-sampled and filtered by the FIR AAF 424 
before being down sampled by the DSB-Q 426. The AAF 424 filters out most or all high 
frequency signals above the receive band. The down sampled signal at the output of the 
DSB-Q 426 may be denoted by z(m). The down sampled signal z(m) may then be filtered 
by the ECF 428 where the output of ECF 428 may be expressed as 

N 

zl{m) = (m - A: + l)wk = z/w, 

where Zm is the data vector 
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Zm [zim\ z(m ~ 1), z(m - 2), ...z(m-N+ l)f, 
w is the ECF coefficient vector, 

A receive signal from Analog Front End may be filtered by Receive Filter 430. 
The output zl(m) of ECF 428 may then be subtracted from the incoming receive signal 
y(m) to generate error signal e(m), e.g., 

e(m) = y(m) - zlOn), 

The samples of the error signal may be generated at the rate i?^. This error signal 
is then used to adaptively train the coefficients of ECF 428 as described in the following 
section. 

According to an embodiment of the present invention, an ECF input z(m) may not 
be the transmit signal but rather a down sampled version. Thus, to derive the LMS 
update rule for the RCDR-EC structure of the present invention, the error signal e(m) 
may be expressed in terms of the down sampled signal z(m\ the receive signal j(m), and 
ECF coefficients as follows: 

^(m)=y(m)-z/w, 

where y(m) is the received signal, Xn is the data vector, and w is the ECF 
coefficient vector as previously discussed. 

Once the error signal e(n) has been defined, the LMS update rule may be given as 
follows: 

"Wm+i =yvm'-jue(m)zn. 
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where fx is the LMS step size. Time subscript for vfrn may represent the LMS iteration. 

The above LMS update rule of the present invention may be similar to a LMS 
update except for the fact that the data vector is the down sampled version of the transmit 
signal x{n), for example. Hence, the computational complexity of LMS update for 
RCDR-EC may be similar as that of the conventional LMS. 

According to another embodiment of the present invention, a dual rate echo 
canceller for simultaneous support of multiple annexes of G.992.1 (or ITU G.dmt) is 
provided. Other annexes and variations thereof may also be supported. The present 
invention provides an efficient Reduced Complexity Dual Rate Echo Chancellor (RCDR- 
EC) architecture for ADSL CPE for simultaneous support of a plurality of annexes, such 
as Annex A and Annex B of G.992,1. In addition, low complexity LMS update rules are 
proposed for adaptive training RCDR-EC. 

The RCDR-EC architecture and LMS update rules of the present invention 
provide various features and advantages. The RCDR-EC structure of the present 
invention may be designed for digital echo cancellation for full rate ADSL at the CPE, 
for example. The proposed structure of the present invention may also support a plurality 
of annexes, such as at least both Annex A and Annex B of ITU Standard G.992.1, also 
known as G.dmt. Other standards, including wire-line technology, may be supported for 
added versatility and improved functionality. In addition, RCDR-EC of the present 
invention may operate at the transmit rate, the lower of the two rates, thereby requiring 
less computation per data sample. 

According to one example, ±e present invention provides an echo canceller 
structure for full rate ADSL CPE. Another aspect of the present invention is directed to 
designing the echo canceller flexible enough to simultaneously support, at least, Annex A 
and Annex B of G.992.1 . Another aspect of the present invention is further directed to 
obtaining an efficient rate matching echo canceller implementation for full rate ADSL 
CPE as well as other applications. 
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For example, the ADSL system, as specified in ITU standard G.992.1, may 
require full-duplex transmission over a single wire line thereby necessitating the 
separation of transmit and receive signals at the ADSL transceivers, such as the CPE. 
Typically, a passive analog circuit, e.g., a hybrid circuit, may be used for this purpose. In 
any realistic implementation, the transmit signal often leaks into the receive path 
whenever perfect balance is not achieved in the hybrid circuit. The part of the transmit 
signal that leaks into the receive path is generally referred to as the echo signal. Even 
though a four-to-two wire conversion circuit has worked well for telephony application, 
the echo signal, if not adequately suppressed, may severely degrade the performance of 
full-duplex ADSL systems. Digital echo cancellers are conmionly used at the ADSL 
transceiver to suppress the echo to an adequate level and to recover the performance loss 
due to the echo. 

ADSL systems support asynmietric transmit and receive data rates, where the 
echo canceller (EC) implemented in the system may be required to accommodate the rate 
difference. For ADSL systems, the upstream rate is generally higher than the 
downstream rate. In other words, for ADSL customer premises equipment the receive 
rate is generally higher than the transmit rate. In accordance with specifications 
recommended by the full rate ADSL ITU standard G.992.1, for example, the following 
may apply: 

where M is an integer larger than 1, Rrx and Rtx are the sample rates of receive and 
transmit signals, respectively. The transmit bandwidth requirements may be different for 
region specific Annexes, e.g.. Annex A and Annex B of G.992.1, for example, and the 
value of M may depend on the particular Annex and/or other factors and conditions. 

The present invention provides a detailed implementation of Reduced Complexity 
Dual Rate Echo Canceller (RCDR-EC) structures for simultaneous support of, at least, 
Annex A and Annex B at the ADSL CPE, for example, as well as other applications. In 
addition, LMS update rules may be implemented for adaptive training of the RCDR-EC 
of the present invention. 
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5— I 



For ADSL CPE and other applications, where the receive rate is higher than the 
transmit data rate, the usual approach may involve up sampling the transmit signal to the 
match the receive rate before echo cancellation. As shown in FIG. 2 above, an up- 
sampled transmit signal may be filtered by the ECF to produce an output signal at the 
5 receive rate. 

However, for such implementation, the ECF may operate at the higher receive 
rate and as a consequence the FIR ECF may require a larger number of taps to span a 
fixed length of time. Thus, extra hardware or software complexities may be needed. 
RCDR-EC of the present invention may operate at the lower transmit rate thereby 
10 reducing hardware and/or software complexities. 

FIG. 5 is a block diagram of a CPE-side echo canceller, according to an 
embodiment of the present invention. According to one example, the RCDR-EC may 
include a combination of a delay block 514, an Input Up Sampling Block 516, an Echo 
Q Cancellation Filter 520, an Output Up Sampling Block 522, an Interpolation Filter 526, 

|T 15 an Anti- Aliasing Filter 532, and a Down Sampling Block 534. In addition, annex 
selectors 518, 524 and/or 536 may be included as well. An optional DC offset block 528 
M= may be incorporated. 

12 A transmit signal may be received as an input to Inverse Fast Fourier Transform 

ffl (IFFT) filter 510 for generating a filtered transmit signal. The filtered transmit signal 

y, 20 may be received by a Transmit Filter 512. A delay block 514 may receive the transmit 
signal for generating a delay in the input transmit samples. 

According to an embodiment of the present invention. Input Up Sampling Block 
(USB-I) 516 may up-sample the delayed transmit signal by a factor of 2 or 4 for Annex A 
and Annex B, respectively. The delayed transmit signal may be up-sampled by other 
25 factors as well, which may vary depending on the applications and/or other conditions. 
The up sampling function may be achieved by zero filling. Other operations may be 
implemented for achieving the up sampling function. The output rate of the USB-I 516 
may be the transmit rate Rtx^ 
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Annex selector 518 may include a number of switches (e.g., three switches) for 
selecting a variety of annexes (e.g.. Annex A or B). As shown in FIG. 5, by way of 
example, the switch position A and B denotes the respective Annexes. Other selections 
and/or options may be implemented. 
5 Echo Cancellation Filter (ECF) 520 may include an adaptive FIR filter that 

receives the up-sampled transmit signal x(n) and generates a copy of the echo, as shown 
by xl(n). As both the input x(n) and output jcl(n) of the ECF 520 may be sampled at the 
lower transmit rate {Rtx\ the ECF filter 520 may operate at the transmit rate Rtx. 

Output up sampling block (USB-0) 522 may receive the output xl{n) of the ECF 
10 520 and up-sample xl(n) by a factor of M. For example, the value of M may be 8 for 
SI Annex A and 4 for Annex B, where the selection may be made via annex selector 524. 

CI Other values of M may be assigned for other options. Other variations may be 

'r^i implemented. The up sampling function may be achieved by zero filling or other 

operation. 

4« 15 The Interpolation Filter (IF) 526 may be a fixed FIR filter that interpolates the up- 

sampled signal x2(m). Essentially, the IF 526 may be a low-pass filter that filters out 
most or all M-l replica of the transmit signal generated due to up sampling. The IF 526 
Oj may operate at the receive rate Rrx and generate an output x3(m) at the receive rate 

The output x3(m) of the IF 526 may be subtracted from an incoming receive signal y{m) 
20 to generate a residual echo signal er{m). 

Anti-Aliasing Filter (AAF) 532 may include a fixed low pass FIR filter that 
receives the residual echo signal er{m) and filters out most or all the high frequency 
signals above the transmit band before feeding the residual echo signal erim) to the down- 
sampling block (DSB) 534. AAF 532 may operate at the receive rate Rr^^ 
25 Down Sampling Block (DSB) 534 may down sample the output Crjim) of the AAF 

532 by a factor M and generate an error signal e(n) at the transmit rate Rtx- For example, 
the value of M may be 8 for Annex A and 4 for Annex B, where the selection may be 
made via annex selector 536. Other values for M may be implemented for other annexes. 
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For RCDR-EC of the present invention, the up sampling by a factor M (e.g., as 
performed by USB-O 522) may be performed after the ECF 520 as opposed to the 
conventional approach where up sampling is performed before the EEC 520. 

According to another embodiment of the present invention, the basic signal flow 
5 within RCDR-EC may be accomplished as follows: The transmit signal x(ri) may be first 
delayed and then up-sampled. In the absence of the delay block 514, the ECF 520 may 
contain a leading sequence of zero coefficients, which may result in an inefficient use of 
the filter coefficients. The leading zero coefficients may attempt to model the delay of 
the echo signal with respect to the transmit signal caused by various elements in the 
10 transmit path and the twisted pair line. The delay block 514 may be added to compensate 
H for this delay of the echo signal and, hence, alleviate the burden for the ECF 520 to 

Q model the delay. According to an example of the present invention, the delay block 512 

^ is not an absolute requirement but may be added for improved efficient use of the ECF 

RJ coefficients. 

>: 15 For example, the highest sub-carrier for Annex A and B may be approximately 

^ around 138 and 276 kHz, respectively. As a result, even though the width of the transmit 

RJ bands may be identical (or similar). Annex B transmit signal may require approximately 

twice the sampling rate for a passband sampling structure. The input up sampling block 
O (USB-I) 516 along with the annex selector 518 may be used to accommodate various 

20 sampling requirements. The output of the IFFT 510 may be sampled at a fixed rate, e.g., 
276 ks/s (kilo samples/sec), for example. After the USB-I 516, the transmit sample rates 
(Rtx) may be approximately equal to 552 ks/s and L104 Ms/s for Annex A and B, 
respectively. The up sampling function may be performed by inserting zeros (e.g., 1 for 
Annex A and 3 for Annex B) in between consecutive samples of the input signals. Other 
25 operations may be used to perform the up sampling function. 

The up-sampled transmit signal may then be filtered by a N-tap FIR ECF 520, 
where x(n) may represent the up-sampled transmit signal and {wi, W2, ws, ... wm}, may 
represent N adaptable ECF coefficients. The output of the ECF 520 may then be written 
as: 
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/V 

xl{ri) = {n'k-\- \)wk = xjw, 



where Xn is the data vector at time n represented by: 



x„ [x(nX x(n - 1), x(n -2\ ...x{n-N + 



w is the coefficient vector represented by: 



W [Wi, W2, ... y^Nf'y 



and T denotes the transpose operation. 

The USB-O 522 may up-sample xl{n) by a factor M. In other words, for every 
input sample, the USB-O 522 may generate M output samples. As mentioned before, the 
value of M may depend on (or be related to) the Annex type and/or other factors, e.g., 



8, for Annex A 
4, for Annex B. 

The up sampling function may be achieved by a zero filling operation (or other 
operation) and the output of the USB-O 522 may be expressed as: 

\x\{n\ ifm = Mn 
x2(m) = < 

[0 otherwise 

In other words, x2(m) may include samples of xl(n) with M-l zeros inserted in 
between. For example, samples of x2{m) may be expressed as follows: 



•0,xl(n-l),0,0,---,0,xl(n),0,- 

M -1 zeros 
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Different time indices (e.g., n and m) may be used to represent samples with 
different rates. For Annex A, n may denote a sampling rate of approximately 552 ks/s 
and for Annex B approximately 1.104 Ms/s. For both Annexes, m may denote a 
sampling rate of approximately 4.416 Ms/s. 

The output of A(/-tap FIR IF 526 may be represented by: 

x3(m)= 2x2(m-fe+l)/^, 

where 

are NflF 526 coefficients. 

The output of IF 526 may be subtracted from the incoming receive signal y(m) to 
generate the residual echo signal, e.g., 

er{m) = y{m) - x3(m). 

The residual echo signal er(m) may be filtered by the FIR AAF 532 to generate a 
filtered version ejf, which may be expressed as: 

where 

{gugl. ^3, ... gNg} 

are Ng AAF 532 coefficients. 

The output of the AAF 532 may be down sampled using DSB 534 by a factor of 
M to generate the error signal In this case, the error signal e(n) may be expressed 
as: 

e(n) = erf(m), where m = Mn. 

In other words, e{n) may include decimated samples of e^n) where the samples 
of e(n) may be written as 

erfiMin - 1)), erfiMnX erf(M{n + 1)), ... 
The value of M (e.g., 8 or 4) may be selected by the Annex selector 536. 
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According to another embodiment of the present invention, the error signal may 
be used to adaptively train the coefficients of ECF 520. 

To describe the LMS update rules for RCDR-EC structure of the present 
invention, the error signal e(n) may be expressed in terms of the transmit signal x(n), the 
receive signal y{m), and IF 526 and AAF 532 filter coefficients. The error signal e{n) 
may be written as: 

e(n) = yr{n)-h^XnW, 

where yr{n) may be the received signal component after AAF 530 and DSB 532, w is the 
ECF 520 coefficient vector, 



x(n) x{n - 1) 
x{n - 1) x{n - 2) 

x(n - L + 1) x(n - L) 



x(n-N + l) 
x(n-N) 

x(n-L-N + 2) 



is the L X M input data matrix, and 

h [huh2,h3, hif, 

is the filter coefficients vector that represents a combined effect of IF 526 and AAF 532 
along with effects of USB-O 522 and DSB 534. The {Nf + A^^ - 1) coefficients of the 
combined filter response of IF 526 and AAF 532 may be represented by bi. For example, 
hi may be obtained by convolving// and gi as follows: 



bi=^f,gi-k^u for/=l,2,...,iV/+N^-l. 
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From the combined filter coefficients b/, the elements of h may be easily obtained. 
The effect of USB-O 522 and DSB 534, in one respect, is to select a phase out of possible 
Af phases. In other words 

hi = bM{i-i)+p. for /= 1, 2, ...,L, 

where j9 % {0,1,2, Af - 1} denotes the selected phase and L may include the smallest 
integer not less than (Nf-^Ng- l)/M, Though the exact value of p may not significantly 
impact the performance, an a priori knowledge of p is useful in calculating h. 
The LMS update rule may be defined as follows: 

w„+i =w„-//e(n)X/h, 

where // is the LMS step size. The time subscript for Wn may represent the LMS iteration. 

At a first glance, it may appear that the proposed LMS update requires (L + l)Ar 
multiplications per iteration. The factor L may be due to the vector-matrix product X^^h. 
However, due to the shift structure of X^, where X^j is a Hankel matrix, for example, per 
iteration, N multiplications may be required. Hence, the over all computational 
complexity of the LMS for RCDR-EC of the present invention may be 2N multiplications 
per iteration. The multiplication required may be approximately twice that of usual LMS 
which in comparison requires only N multiplications. Since, there is a saving of filter 
coefficients by a factor M, the computational complexity of the proposed LMS for 
RCDR-EC of the present invention is lower for at least both values of M = 4 or 8. 

A simplified LMS may be used to reduce computational complexities by using an 
approximation on h. 

As discussed above, the IF 526 and/or AAF 532 may include low-pass filters. 
Generally speaking, the impulse responses may be similar to the sine function. In other 
words, the impulse response may exhibit a dominant peak and diminish in magnitude 
away from the peak. Moreover, the convolution of the two such filters may also exhibit a 



-37- 



PATENT 

Attorney Docket No. 56162.000322 



similar behavior. Other types of filters may be used for IF and/or AAF, in accordance 
with the present invention. 

For example, both IF and AAF may include low-pass filters for filtering out the 
high frequency signals above the transmit band. According to one example, the same 
5 (similar or related) coefficients may be used for both IF and AAF, e.g.,/^ = gk for all k. 
For an Annex B application, a 40-tap low pass filter design may be used for at least one 
of IF and AAF, for example. The impulse response of the 40-tap filter along with the 
time-domain convolution of the two filters are illustrated in FIG.s 6 and 7. FIG.s 6 and 7 
show that the resulting filter after convolution exhibit the properties discussed above. In 
10 particular, FIG. 6 illustrates an impulse response of IF and AAF. FIG. 7 illustrates a 
convolution of IF and AAF. 

■y According to an embodiment of the present invention, the largest magnitude 

d 

ru element of h may be replaced with its sign. Also, the rest of the elements may be set to 

S'j zero. Replacing the largest element with its sign may involve the scaling of h where such 

15 scaling may be absorbed in ft. Accordingly, there may be little or no approximation 

g: involved. After such scaling, the remaining elements having a magnitude less than one 

iZ may be approximated to zero. With this approximation the coefficient vector h becomes 
a vector of zeros except one non-zero element (± 1). Without loss of generality, the 

P'J largest magnitude element of h may be assumed to be positive. The matrix-vector 

20 product may be reduced to a delayed version of the input vector, e.g.. 



An h = 



x(n - d) 
x{j%-d-i) 

x{n-d-N + l) 



25 where h is the approximate h and d is the index of the largest element of h. 
Using this approximation, the LMS update may simplify to, 
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w„+i = Wn-//e(«)x«-c/, 
where x„ is the data vector defined before as 

5 

Xn = [x(n), x(n - 1), x{n - 2), . . . jc(n - + 1)]^. 

Hence, the simplified LMS update for RCDR-EC is similar to the normal LMS 
update except for the delay. In other words, the computational complexity of the 
10 simplified LMS for RCDR-EC of the present invention is similar to that of usual LMS. 

Since both IF and AAF are part of the system design, the approximation does not 
restrict the use of the simplified LMS. Thus, the IF and AAF may be designed such that 
the resulting combined filter has at least a single dominant coefficient. Factors that may 
severely degrade the convergence and performance of LMS may include multiple 
15 dominant peaks in the combined filter response and inexact setting of the delay d. 

The USB-0 522 with the zero filling structure may be combined with the IF 526 
and efficiently implemented using a poly-phase structure, for example. The same may 
apply for the DSB 534 and AAF 532. 

Differences between conventional echo canceller and RCDR-EC of the present 
20 invention may include various factors, such as structural differences, filter lengths, 
computation for LMS training, stead-state computational complexities and noise floor 
levels, for example. 

As for structural differences, a conventional echo canceller may include an up 
sampling block followed by the adaptive FIR ECF that operates at the higher receive rate 
25 Rrx> In the RCDR-EC approach of the present invention, an adaptive FIR ECF may 
operate at the lower transmit rate Rtx^ In addition, a USB and a IF may be used to 
generate the echo samples at the receive rate and a DSB and an AAF may be used to 
calculate the error signal at the transmit rate. 
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As for filter lengths, to capture the same (or similar) fixed amount of time span 
(e.g., echo tail), the RCDR-EC of the present invention may include N filter taps as 
compared to the MN taps required by the conventional echo canceller approach. 
However, in addition to the ECF, RCDR-EC may implement at least two more filters, 
5 e.g., IF and AAF, for example. In general, the ECF is longer than the combined length of 
IF and AAF and hence, there is a significant reduction in number of filter taps. For 
example, for a typical value ofN = 256 and M = 8, the conventional filter may require 
approximately 2048 taps. On the other hand, even with two 64 taps IF and AAF filters, 
the number of taps for RCDR-EC of the present invention is approximately 384, 
10 substantially less than that required in o±er systems. 

As for LMS training computation, the computational complexity may be similar 
p for the conventional LMS and simplified LMS for RCDR-EC of the present invention. 

rt'J For example, as the conventional ECF may be M times longer than that of the RCDR-EC 

D of the present invention, the number of multiplications required for LMS update may be 

in I 

15 M times higher for a conventional echo canceller. Due to the reduction in the filter taps, 
even the exact LMS for RCDR-EC may have M/2 times lower complexity then the 

H conventional approach. 

"-^1^ As for steady-state computational complexity, in steady-state using efficient poly- 

^ phase implementation, for example, the total number of multiplications per second 

20 associated with RCDR-EC of the present invention may be represented by 

In contrast, even with efficient poly-phase implementation, the number of 
25 multiplications per second required by the conventional echo canceller is 
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In a typical implementation, the total filter length of IF and AAF is less than the 
filter length of ECF, e.g., Nf + Ng<N. As a result, RCDR-EC of the present invention 
may reduce the steady-state computation by at least a factor M. 

As for noise floor level and overall performance, in a RCDR-EC system, the AAF 
filters out signals above a transmit band before generating the error signal e{n). Thus, the 
error signal reduces the contribution from the wide band receive signal which acts as 
noise for echo cancellers. In a conventional echo canceller, most or all the out of band 
receive signal works to contribute to the error signal. According to the present invention, 
this reduced noise floor for RCDR-EC results in reduced convergence time and improved 
overall performance of IMS. Other comparisons as well as advantages may be 
recognized by the present invention and variations thereof. 

According to another embodiment of the present invention, a weighted error echo 
canceller for transceivers with unequal bandwidths is provided. The present invention 
enables an adaptive echo cancellation filter to perform at a lower sample rate. 

Echo cancellers (EC) may be conmionly used in full duplex conmiunication 
systems to suppress the echo of the transmit signal entering the receive path. When the 
transmit and receive path have different bandwidths (hence, different sample rates), the 
echo canceller may need to bridge the rate difference. With conventional designs, the 
echo canceller filter operates at the higher sample rate. According to the present 
invention, an efficient echo cancellation scheme, referred to as Weighted Vector Error 
Echo Canceller (WEVE-EC) is proposed for various applications, such as applications 
where the receive path has a higher sampling rate. The WEVE-EC architecture may be 
implemented along with an adaptive algorithm, which may be based on Least Mean 
Square (LMS) update rules. Various other adaptive algorithms may be used to train the 
WEVE-EC of the present invention. 

According to an embodiment of the present invention, an error vector, rather than 
a scalar error signal, may be implemented for updating the echo canceller. In a 
conventional approach, essentially all the sampling phases of the error signal are treated 
as consecutive time samples of a scalar signal at the receive rate. In the proposed 
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Structure of the present invention, most or all the sampling phases of the error signal may 
be stacked in a vector and treated as a vector sampled at the transmit rate. 

An Error Weighting Multi-input-multi-output Filter (EWMF) of the present 
invention provides a flexible weighting scheme on most or all the sampling phases of the 

5 error signal. The Multi-Input Multi-Output (MIMO) filter structure is more general than 
the conventional Linear Time-Invariant (LTI) Finite Impulse Response (FIR) filters. In 
particular, the LTI FIR structure may be considered a special case of the proposed 
EWMF. Hence, the proposed structure of the present invention provides greater variety 
in selecting a training method for echo cancellers. 

10 According to another embodiment of the present invention, the combination of an 

adaptive filter and an interpolation filter enables the echo canceller to operate at the lower 
transmit sample rate. This is in contrast to the conventional echo canceller approach 
where the echo canceller operates at the higher receive sampling rate. In addition, 
training and updating functionality of the echo canceller may be conducted at the same 

15 (or similar) low sample rate. 

The architecture of the present invention reduces the number of filter coefficients 
without loss of performance. As a result, hardware requirement and computational 
complexities are significantly reduced and an increase in convergence speed of training 
algorithms may be achieved. Other advantages may be observed from the various 

20 embodiments of the present invention. 

The proposed architecture may be suitable for steady-state tracking of the echo- 
canceller. The proposed architecture may also be easily modified for other situations, 
such as the reverse situation where the receive bandwidth is smaller than the transmit 
bandwidth. 

25 The present invention provides an effective echo cancellation method that takes 

advantage of the situation where the transmit signal, and thus the echo signal generated 
from it, has a smaller bandwidth than that of the receive signal. The architecture of the 
present invention may include various components and function blocks, such as an 
Adaptive Echo Cancellation Filter (AECF) 810 which may operate at the lower transmit 
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sample rate, an Interpolation Filter Bank (IFB) 820, and an Error Weighting MIMO Filter 
(EWMF) 830. Other functions may be implemented in the system of the present 
invention. 

A corresponding LMS algorithm of the present invention provides training, 
5 tracking and/or other purposes. A filter bank design provides an efficient implementation 
for enabling some or all functional blocks to operate at the lower sample rate. Compared 
with the conventional echo cancellation filter which operates at the higher sample rate, 
the structure of the present invention offers improved performance, while significantly 
reducing the required number of taps, the sample rate, as well as yielding improved 
10 training behavior and less power consumption. Moreover, the proposed algorithm of the 
. , present invention also yields better convergence and stability properties during training 
C and other phases. 

pj FIG. 8 illustrates a basic signal flow and the functional blocks of the WEVE-EC, 

S'; according to an embodiment of the present invention. The Weighted Vector Error Echo 

a s i; 

N=15 Canceller (WEVE-EC) of the present invention may include various functional blocks, 
T' such as an Adaptive Echo Canceller Filter (AECF) 810; an Interpolation Filter Block 
^: (IFB) 820; an Error Weighting Multi-input-multi-output filter (EWMF) 830; a Reference 
U Signal Interpolation Filter (RSIF) 840; and a Vectorization Unit (VU) 850. Otiier 
Si functional blocks may also be implemented, in accordance with the present invention. 
N^20 The AECF 810 may include a FIR filter with Nw coefficients and an adaptive 

component of the WEVE-EC, The AECF 810 may operate at the lower transmit sample 
rate. The A^^ coefficients of the AECF may be stacked into a vector w, where 

25 

The coefficient vector w may be adapted during a training period and may be 
updated during the steady state. The AECF 810 may receive the transmit signal x„ and 
generate a copy of the echo signal that is fed to the IFB 820. 
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The JFB 820 {fu ... , /m} may include a fixed filter bank for generating M 
branches of signals firom the output of the AECF 810 to match a receive signal vector. 
The branches of signals may be organized into the vector y„ = |>„' , . . . , yj^f at each time 
instance or other time interval. For example, an integer M may be used where the 
architecture of the present invention may be easily modified with a scheme, such as a 
mixed upsample-downsample scheme, to accommodate various rational values of M. 
The output of the IFB 820 may include a vector signal which may be obtained by 
filtering the output of the AECF 810. This signal may be used to generate the residual 
echo Tn. 

The weighted error vector represents the training error of the adaptive echo 
canceller filter. The EWMF 830 H allows a flexible weighting on the error vectors and 
thus offers a vast set of options for advanced training of the WEVE-EC of the present 
invention. 

The RSIF 840 may include a fixed filter for generating a reference signal for 
trainmg and/or other purposes. The filtering operation may involve a combination of the 
IFB 820 and the EWMF 830, for example. 

The ti:ansmit signal x„ may occupy a different (e.g., smaller) bandwidtii than tiiat 
of tiie receive signal. The filtered tiransmit signal q) may be obtained by filtering the 
ti^smit signal x„. The signal x„ may then be used to generate the reference signal s). 

The reference signal s) may be obtained by filtering q). The reference signal may 
be used for ±e LMS algorithm during the tiraining of die AECF 810 with coefficient 
vector w. 

The receive signal vector d„ may be represented as 

This receive signal vector may be obtained by stacking the M sampling phases of 
tile receive signal into a M length vector. The receive patii may provide a sample rate M 
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times higher than the transmit path. Hence, the vectorized receive signal d„ has a 
sampling rate equal (or similar) to the transmit sampling rate. 

The VU 850 may include sub blocks, such as Delay Units z(k) and Down 
Sampling blocks. Other sub blocks may also be incorporated by VU 850. The Delay 
Units sub block as designated by z(k) may delay a signal by K samples. The Down 
Sampling blocks as designated by M and a down arrow may be used to drop M - 1 
intermediate samples from the input signal. In other words, Down Sampling block 
decimates the input samples by a factor M. 

The residual echo vector r„ = d„ - y„ may be obtained by subtractmg the output of 
the IFB 820 from the receive signal. During the training period, the residual echo vector 
r„ may be used to generate the weighted error vector e„ for the LMS update or other 
purpose. The weighted error vector e„ may include the filtered residual echo vectors and 
may be used for training purposes, e.g., to train the AECF 810. Once the WAVE-EC is 
properly trained, e„ may contain a small portion of the echo signal. 

In the steady state, operational function blocks may be limited to the AECF 810, 
the IFB 820, and the VU 850. The RSIF 840 and the EWMF 830 may be implemented 
for training/tracking purpose. The signals for steady state operation may include the 
echo-canceller output y/ and the echo cancelled signal r/. The signals Sk and ek may be 
generated for the adaptive LMS. 

Various echo canceller parameters for the architecture of the present invention as 
well as the LMS algorithm may be implemented as described in detail below. 

Liteger M may include the ratio of the receive sample rate to the transmit sample 
rate. M > 1 implies that the transmit sample rate is lower than the receive sample rate. 
According to the present invention, the IFB 820 may be used on the transmit signal and a 
branching operation on the receive signal to bridge the transmit path and the receive path. 

Adaptive Echo Canceller Filter 810 length may represent the number of taps in 
filter 810. As the architecture of the present invention may operate at a lower sample 
rate, the length is substantially lower than that of an echo cancellation filter operating at a 
higher sample rate. 
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EWMF coefficients may be grouped into a matrix, such as {hm^nl m = 1, M, n 
= 1, Nh}> The choice of the coefficients may depend on the nature of the transmit 
signal, receive signal, and/or the echo signal, as well as, the training algorithm used to 
train the echo canceller. Other factors may also be involved in the choice of coefficients. 

The dimension of the EWMF matrix may be represented as M x Nh^ For LMS 
algorithms, Nh may be a multiple of M, e.g., K = Nh/M is an integer. K may represent the 
number of residual echo vectors that are used to generate the weighted error vector. A 
larger EWMF matrix permits a more flexible weighting scheme on the error signal at the 
cost of higher complexity. 

Interpolation Filter Bank (IFB) coefficients may be represented by {fn^^\ ik = 1, . . ., 
M, n = 0, . . . 4 } : A dedicated filter of length Nf may exist for each of the M signal 
branches. These coefficients may be chosen such that a desired interpolation filtering 
performance is achieved with a given filter length. 

Interpolation Filter Bank length may be represented by Nf. The length may be 
decided by a desired performance of the up sampling block, the computational cost 
and/or other factors. The performance as well as the complexity may increase with the 
length of the filter. 

LMS update step size sequence may be represented by Generally, the larger 
the step size, the faster the training will converge. However, there may be a limit on the 
magnitude of the step size. For example, a large enough value may cause the training 
algorithm to diverge. On the other hand, the smaller the value of fin. the smaller the error 
floor. For example, the step size sequence may be chosen as constant or as a decreasing 
set of numbers such that the algorithm quickly achieves convergence with large step sizes 
at the early stages of the training procedure and obtains low error floor with smaller step 
sizes later in the training procedure. This is conmionly known as a gear shifting method. 

The training procedure of the present invention may involve adapting the 
coefficients of the echo canceller filter based on up^sampled reference sequences and/or 
weighted error signal. Other factors may be considered. The adaptive echo canceller 
may be matched to the echo path such that a significant component of the echo 
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component in the receive signal may be cancelled. The effectiveness of the echo 
cancellation may be measured by the energy of the weighted error signal. The tracking 
procedure may use the same (similar or related) update mechanism as that of the training 
procedure or a variation thereof. 
5 The training method is not limited to the LMS algorithm as discussed above, but 

other algorithms may be used instead of LMS, in accordance with the present invention. 
The scalar form of an LMS update rule may be described as follows: 

10 

M. where Wn is the adaptive entity for the LMS algorithm, e« is an error signal, n is 

i;; the time index, and /Un is the LMS step size. The update rule for the proposed LMS 

jli algorithm is shown below in step 918. Compared with the scalar form, the proposed 

update rule of the present invention provides a more flexible weighting scheme for the 
15 error signal. 

B In the proposed LMS algorithm, each update may involve the steps of FIG. 9. At 

step 910, a vector corresponding to a reference filter may be formed. At step 912, a 
H matrix may be formed. At step 914, a weighted error vector may be formed. At step 916, 

p an adaptive echo canceller filter may be updated. At step 918, the weighted error signal 

^ 20 energy may be calculated. The steps of FIG. 9 will be described in detail below. 

At step 910, vectors 8;^^ ... , corresponding to reference filter may be formed 

where 

s/=[/„ S^n-l T s\.M^+lf. 

25 

At step 912, the matrix may be formed as follows: 

S/i = [ Sn T Srt^ ]. 
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At step 914, weighted error vector may be formed as follows: 

e„=[e„' e„' t e/f. 
At step 916, the adaptive echo canceller filter may be updated as follows: 

Wn+l =Wn+j«„Sne„. 

At step 918, the weighted error signal energy e„\ may be calculated. At step 
920, it may be determined whether the energy is less than a predetermined value. If so, 
the process may be terminated at step 922. If the energy level is equal to or above a 
predetermined value, then n = r + Land a feedback loop to step 910 may be estabhshed. 

The present invention also provides advantages with respect to computational 
complexity. For example, it may take N,^ + MNf + MNh multiplications to generate the 
weighted error vector Cn, MNf + MNh multiplications to form the matrix S„, and M^Ny, 
multiplications to execute the update equation for each update of the AECF. In total, (1 + 
M'^)Ny„ + IMiNf+Nu) multiplications may be required for each update of the AECF. 

Additional details of the LMS algorithm of the present invention will now be 
described. The output of the echo canceller interpolation filter bank (e.g., IFB 820) 

may be written as 

y„ = FX„w„ 

where F represents the echo canceller Interpolation Filter Bank, e.g.. 
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f(0 f(l) 

Jo J\ JNf-l 

f(2) r{l) nil) 

P _ ^0 J\ JNf-l 

f(M) MM) jriM) 
Jo Jl JNf-\ 

and 

'n-Nf-Nw+2_ 

is the transmit signal matrix. 

In one example, the filter bank F may implement a zero-insertion plus low-pass 
filtering operation. In such cases, F may be written as 





fo 




f{Nf-\)M 


F = 




fM+l 


f{Nf~X)M+\ 




_/m-i 


flM-\ 





where [fo, . . . jNjM-\f may be a iV/M-length low-pass filter. 

As the receive signal may occupy a larger bandwidth than that of the transmit 
signal, the residual echo vector may have a significant out-of-band noise. This out-of- 
band noise has detrimental effect on the adaptation process and should be removed. 
According to an embodiment of the present invention, the out-of-band noise may be 
removed or minimized by passing the residual echo vector through the EWMF, or other 
filter. As discussed above, the residual echo signal may be represented by 
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r« = d„ - y„. 

Using the residual echo signal, the weighted error vector may be written as 















= H 





In the formulation, K = Nr/M residual echo vectors may be involved in the 
filtering operation to generate a filtered error vector. By using the matrix 



5 

to represent the weighting operation, the possibility of using different filters for each 
sampling phase is preserved. For the special case where a single K tap FIR filter is used 
for all sampling phase, the H may become a Toeplitz matrix with the first row 

M-\ zeros M-1 zeros Af-1 zeros 

where the jRTnon zeros elements {hu ,,,,hK} maybe defined as the FIR filter coefficients. 
A Toeplitz matrix generally refers to a matrix whose entries are constant along each 
diagonal. 

After defining the necessary signals, the corresponding Least Mean Square (LMS) 
algorithm for the coefficients of the echo cancellation filter may be derived in accordance 
with the present invention. For example, the instant gradient for the LMS training 
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algorithm may be obtained by taking the derivative of e„^e„ with respect to w. The 
update rule may then be represented as: 



2 dw„ 



= w„ + M« [Xjr I T I XVk.iF^] • • e„. 

The matrix multiplication in the update rale may represent the operations 
performed in the reference signal interpolation filter (e.g., 840) in FIG. 8. 



[X„^F^|...|XVk^,F^]-H 



M 



M 



M 



where the filtered output s/ may be arranged into a matrix. The final form of the update 
rule may be represented as 
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According to another embodiment of the present invention, a minimum mean 
square error linearly constrained fast algorithm for adaptive training of a Time Domain 
Equalizer (MLC-TEQ) is provided. A fast adaptive algorithm of the present invention 
may be used to obtain Finite Impulse Response (FIR) filter coefficients for Time domain 
Equalizer (TEQ) used in Discrete Multitone (DMT) based applications, such as ADSL, 
for example. The TEQ coefficients obtained by the algorithm of the present invention 
shortens the overall effective discrete time channel impulse response length within a 
given target length (e.g., symbol prefix length for DMT application). Advantages of the 
proposed data aided adaptive algorithm may include providing the TEQ filter coefficients 
with near-optimal performance; having low computational requirements, having fast 
convergence, and exhibiting attractive stability properties. Other advantages may also be 
realized by the present invention and variations thereof 

Other features of the algorithm of the present invention may include 
implementing a Target Impulse Response (TIR) vector, e.g., the combined response of 
the transmission channel and TEQ filter, as a search vector, e.g., the vector that needs to 
be calculated. 

The algorithm of the present invention may constrain a selected element of the 
TIR vector to be approximately equal to a constant where the selected element may be 
any element of the TIR vector. The error may be computed by using the difference 
between the TEQ output and TIR filter output. TEQ filter coefficients may be adaptively 
computed by minimizing the Mean Square Error (MSE) criterion. TEQ filter obtained 
through the algorithm of the present invention may provide minimal energy Inter Symbol 
Interference (ISI) and Inter Channel Interference (ICI). 

The computation requirement of the algorithm of the present invention may be 
linear in equalizer length, e.g., the number of filter taps of the equalizer. As a result, the 
algorithm is suitable for practical implementation in various applications. The algorithm 
of the present invention may adaptively calculate an equalization delay using the training 
and received sequences thereby providing efficient use of FIR TEQ filter coefficients. 
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The algorithm of the present invention may incorporate input data conditioning through 
use of dummy signals thereby providing enhanced stability and improved convergence. 

The present invention provides a fast, stable and high performance adaptive 
algorithm for computation of near-optimal FIR TEQ coefficients. The resulting TEQ 
filter shortens the overall channel response such that residual energy of the overall 
channel outside the target length is minimized, hi applications, such as DMT, TEQ 
filters obtained through the algorithm of the present invention provides minimum energy 
Inter Symbol Interference (LSI) and Inter Channel (or Carrier) Interference (ICI) and as a 
result yields a higher bit rate. Odier advantages may also be recognized. 

In Discrete Multitone (DMT) systems, it may be desirable for the discrete time 
channel spread to be less than a prefix length to avoid ISI and ICI. In some applications, 
the spread of the channel may be greater than the prefix length. Hence, preprocessing of 
data may be implemented to shorten the effective channel seen by the DMT receiver. An 
approach, e.g.. Time domain Equalization (TEQ), may be implemented to apply (sample 
or fractionally spaced) linear filtering to the sampled receive data to shorten the effective 
channel spread. The present invention provides an approach to designing sample spaced 
and fractionally spaced (e.g., with factor 2 or other factor) Finite Impulse Response (FIR) 
TEQ filters. 

FIG. 10a outlines a basic setup for a TEQ system, according to an embodiment of 
the present invention. Here, Xk represents the ^-th sample of the transmitted signal with 
variance x- Block 1010 with up-arrow and M corresponds to an up-sampling fiinction 
with factor M, which may involve inserting M - 1 zeros between the samples of Xk. For 
example, for M = 2 case, 4 sequence may be given by 

0, X.2, 0, X.l, 0, Xo, 0, Xi, 0, X2, 0, X3, 0, <r- 

For example, M = 1 may refer to a sample spaced TEQ, as represented in Fig. 
10b, and M = 2 may refer to a fractionally spaced TEQ. 
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The combined effects of the transmit filter shaping, the receiver filter, and the 
distortion caused by the transmission channel may be modeled by a Linear Time 
Invariant (LTI) system, as illustrated by 1012, with impulse response denoted by hk. The 
combined effects of the receiver noise and the interferences may be modeled by the 
additive disturbance v*. 

In FIG. 10a, block 1014 represents a TEQ filter having Nteq filter coefficients. 
Nteq TEQ coefficients may be denoted by {w,. W2, w^-eq} and TEQ coefficient 
vector may be denoted by w, e.g.. 



An output of the TEQ filter 1014 may be down-sampled where the down 
sampling process is represented by a block 1016 with down-arrow and M. The down- 
sampling operation may involve selecting 1 out of M samples of oj. For example, for M 
= 2 case, Zk sequence may include 



The time domain equaUzation may convert die channel impulse response hk into 
an impulse response with Nt nonzero consecutive entiies. Therefore, if tiie nonzero 
entries of the resulting channel are modeled with the sequence {Z>*; it % {0 4- ATj- - 1 } } , the 
target channel may be modeled as 



where d is the effective delay corresponding to the starting location of the non-zero 
segment of the impulse response. 



w = [ Wi ^ Wnteq f 



4-0.2, Oo, 02, O4, <— 




n-d 



d<n<d + Nj.-l 
otherwise. 
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In DMT based DSL systems, such as G.dmt and G.lite, standard generally 
dedicates certain training sequences for the computation of TEQ coefficients. These 
training sequences may be periodic DMT symbols without prefix where the period is one 
symbol long, for example. According to one example, {pn, n % {0, <^ , Nr - I}} may 
5 represent the DMT symbol of length Nr, where the transmitted training sequence r„ 
represents the repeated version of pn, e.g., 

r„ = pm and m = mod(n, Nr) 

10 where fnod(n, Nr) is the remainder of n/NR and hence, m is an integer between 0 

P and Nr- I. 

y Therefore, may represent a periodic extension of the sequence p„, e.g., r„ may 

p be obtained through periodic repetition of p„ where the period is Nr. The periodic 
extension may be expressed as 

* 15 

y= PNR-2, PNR-U PO, PU P2, <r- , PnR-U PO, Pi, <r- 

An adaptive TEQ problem may involve developing an algorithm that processes 
training signal and corresponding received signal y„ to obtain TEQ filter coefficients 
20 without any a priori knowledge of channel and noise statistics. 

For DMT applications, an optimal approach may involve choosing {wk} such that 
the number of bits loaded per symbol may be maximized. However, this approach may 
be difficult to implement and impractical for adaptive implementations. Other alternative 
schemes either have high computational complexity and are therefore not suitable for 
25 adaptive implementations or have lower complexity but low performance and instability 
in the convergence. The present invention provides an adaptive approach that yields 
enhanced (or near-optimal) performance. 

MLC-FAST-TEQ algorithm of the present invention uses training data sequence 
{r;j}, and the corresponding received data sequence {yn}. Exemplary special cases for the 
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algorithm may include a Sample Spaced Case where M = 1 and a Fractionally Spaced 
Case with factor 2 where Af = 2. In the Fractionally Spaced Case, {yen} and {yon} may 
represent even and odd phases of the received sequence respectively, e.g., 

yen = yin 
yon = yim-Xy 

where yn is the over-sampled received sequence. 

In this section, various algorithm parameters, which may include predetermined 
constants independent of the data, may be used in the algorithm of the present invention. 
Their values may be adjusted to achieve different levels of performance and stability as 
well as other results. 

The algorithm may assume certain conditions about the amount of training 
sequence data required, e.g., the length of the training sequence Mr- For example, the 
training sequence may involve consecutive samples of the received signal {y^; n % {0, ... , 
Ntr - 1 } } for the sample spaced case. In another example, the training sequence may 
involve consecutive samples of the sequences {y,„; n % {0, ... , AT^, - 1 } } and {y^„; n%{0, 
... , A^rr - 1 } } for the fractionally spaced case. 

Parameters used in the algorithm of the present invention include: Nt which 
represents target channel length; Njeq which represents TEQ filter length; Nr which 
represents symbol length; L which represents fixed tap location where I X {0, ... , ATy - 1 } ; 
Ns which represents TEQ algorithm search vector size Ns = Nt + Nteq; ® which 
represents least square cost function forgetting factor; No which represents maximum 
tone number for the dummy signal and N^^ which represents number of update variables 
where 
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•sr:j» 



N,= 



3 Sampled SpacedCase 

4 FractionallySpacedCaseWithFactorl, 



Other parameters may include {Q; A: % {0, ^ , iVz? - 1}} which represents 
frequency domain dummy signal values. There are No such complex values, e.g., Ck may 
be complex numbers. Other parameters may include { } which represents time domain 
dummy signal, which may be determined by the frequency domain dunmiy signal values. 
The following equation shows the relation: 



>n = Co+ X^RealpartoflC^^^"^;,} nx {0, 4- ,iV^- 1}. 



10 

where the term Real part of stands for the real part of a complex number. The time 
domain dummy signal may be defined as the repeated version of ©n, e.g., {©n = ©m }, 
where m = mod{n, Nr), In other words, ©n may represent the periodic extension of the 
signal ®n. 

15 Other parameters may include S which represents the scaling of the received 

dunrniy signal; c^r which represents Internal dynamic range arrangement parameter; and 
ca which represents a priori variance constant used for the initialization of the matrix F. 

The MLC-FAST-TEQ algorithm of the present invention includes at least two 
steps, which may include the adaptive computation of the equalization delay parameter d 

20 and the adaptive computation of the TEQ filter coefficients based on the delay d 
calculated in the first step. 

An estimate of cross-correlation function {q(k); A: % {0, <- , M x A/i? - 1 } } may be 
computed as 



25 qik)= X ^r^yiMn + k-hlNR)} 

1=0 «=0 
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where NoMNr < Ntr. 

Then ® may be defined as the argument maximizing ±e absolute value of the 
cross-correlation function: 

® = argmaxk \q(k)\ 

Therefore, ® represents the peak point location for the absolute cross-correlation 
function. 

An equalization delay d may be selected as a function of ®. An appropriate 
choice of d may be given by: 



d=l/M{ ® + Nteq/2}. 



In other words, 



d = 



TEO 

K + — — Sampled SpacedCase 

1C TEO 

— + ^ FractionallySpacedCase 



Adaptive computation of TEQ coefficients for a sample spaced case includes 
algorithm variables such as the following: 
^ : Column Vector (Size N;^ x L); 
© : Column Vector (Size N^^ x 1.); 
? : Column Vector (Size N^^ x L); 
® : Column Vector (Size N^^ x L); 
^ : Column Vector (Size N^ x 1.); 
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10 



y 



15 



20 



A : Matrix (Size N^^ x Ns), initialized to all zeros; 
D : Matrix (Size N^ x Ns), initialized to all zeros; 
X : Column Vector (Size Ns x 1), initialized to all zeros; 
k : Column Vector (Size Ns x 1), initialized to all zeros; 

F : Matrix (Size N^^ x N^^), initialized to gjInx, where In^ stands for identity matrix 
of size Nj^ X N^; 

t : Column Vector (Size N^^ x 1), initialized to all zeros; 

b : Column Vector (Size Ns x 1), initialized to all zeros; 

m : Colunm Vector (Size Ns x 1); 

(3) : Column Vector (Size N^^ x 1); 

c : Colunm Vector (Size Ns x 1), initialized to all zeros; 

X : Scalar; 

e : Scalar; and 

n : Scalar. 

TEQ coefficient algorithm steps include the following: 



step 1. 
step 2. 



n = 0, 



'n~d-L 



n-d-L-\_ 



Step 3. 



y 



step 4. 
step 5. 
step 6. 
step 7. 



/b^ = ® +AX, 
A = A-k*^\ 
?=:-^(l-k^X), 
F = (DF, 
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step 8. 

step 9. 
step 10. 
step 11. 



step 12. 



step 14. 



t = F 



1 T 

a 



F = F - t(^'^F), 
t = caF-^, 
b = k + At, 



m = 



\NjEQ+L):(Ns-V) 



step 13. ®'^=[bNTBQ 



9NTEQ4-H 



bNsl. 



X = 



I. 

Y 



step 15. 




® +DX, 




step 16. 


k = 






step 17. 


D = 






step 18, 


A. 

X - 






step 19. 


e = r 


n-d-L-1 + ®n-d-I^l ' X , 


step 20. 


c = c + ke, 




step 21. 


n = n + 1, 




step 22. 


if n < Ntr then go to step 2, else stop 



TEQ coefficients may be given by the first Nteq taps of the vector c, e.g., 
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W = Ci:NTEQ 

Adaptive computation of TEQ coefficients for a fractionally spaced case with 
Factor 2 (e.g., M=2) may include variables that are identical (or similar) to that for the 
sampled spaced case, discussed above. The steps for the algorithm may include the 
following: 



step 1. 



step 2. 



step 3. 

step 4. 
step 5. 
step 6. 
step 7. 

step 8. 

step 9. 
step 10. 
step 11. 



n = 0. 



(D = 



yon 



?=(D+AX, 

F=®F, 

f \ 
e 



t = F 



1 ,r 



+ e'Fe 



a 



F = F-t(.^'^F), 
t = c«FT 
b = k + At, 
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step 12. 



m = 



b. 



(NrBQ+L).(Ns-\) 



Step 13. (D'^ = [bNTEQ 



bNTEQ+H bNs]. 



Step 14. 



X = 



X 



(^7-£e+I).(A^ree+jL-2) 

Y 

^(A/7Eg+L):(A/5-l) 



step 15. 


^ = 


® + DX, 




step 16. 


k = 






step 17. 


D = 


D - k-^^. 




step 18. 


X = 


c'^X, 




step 19. 


e = r 


n-d-L-l + ®n-d-L-l " X , 


step 20. 


c = c +ke, 




step 21. 


n = n+ 1, 




step 22. 


if n < Ntr then go to step 2, else stop. 



TEQ coefficients are given by the first Njeq taps of the vector c, e.g. 

W = Ci:NTEQ 
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Further to the adaptive manner of obtaining FIR coefficients in a TEQ filter 
described above, we now turn to a method and system for shortening the channel impulse 
response in DMT systems wherein a model is employed to effectively represent a 
hypothetical delay transmit signal in conjunction with a desired shortened channel 
impulse response filter to arrive at near optimal TEQ coefficients. In contrast to the 
adaptive process described above, where the TEQ setup involved a fractionally spaced 
TEQ scenario in which up-sampling/down-sampling factor M = 2, as shown in Fig. 10a, 
the TEQ filter solution described in this section is based on a sample spaced TEQ setup in 
which M = 1, as shown in Fig. 10b. 

In DMT systems, each received data symbol or packet is commonly prepended 
with a guard sequence called the cyclic prefix that is used to space the symbols so that the 
receiving system recognizes each individual packet of data. Each DMT symbol typically 
consists of N samples (where N is an arbitrary integer), while the cychc prefix is a copy 
of die last V (where v is an integer less than N) samples of the DMT symbol. In order to 
prevent ISI, the value of v must be greater than or equal to the length of the channel 
impulse response, i.e, the channel length, commonly denoted by the variable h. 

Upon receipt of the symbol by the receiving modem, the cyclic prefix is 
discarded, since it includes no information not otherwise included within the symbol. 
The symbol is further processed, and conventional methods of shortening the channel 
impulse response are applied. Because the cyclic prefix does not convey any new 
information about the transmitted signal, efficiency of the system decreases in proportion 
to N/(N4-v). Therefore, in order to maximize system efficiency, either N should be large, 
or V should be small. However, as N increases, memory requirements in the 
communications system are also increased. 

In DMT systems, it is further desirable that the discrete time channel spread is 
less than the prefix length so as to avoid ISI as well as ICI, which is the convolution of 
signals on overlapping carrier channels. In applied systems, the spread of the channel is 
typically significantly greater than the prefix length. Therefore, preprocessing of data is 
required to shorten the effective channel seen by the DMT receiver, thereby reducing the 
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occurrence of ISI and ICL The most common approach to this preprocessing function is 
to apply Time Domain Equalization (TEQ) to the sampled receive data to shorten the 
effective channel spread. TEQ refers to a process for filtering the incoming signal in 
such a way as to combine the data signals from various channels and reduce the cyclic 
5 prefix of the resulting signal to a uniform size that is as small as possible. 

FIG. 10b illustrates the basic scenario involving the above described TEQ 
problem. In particular, Xk (where /: is an integer) 1000 denotes the transmitted signal that 
is propagated along a channel and ultimately received by a modem. The transmitted 
signal includes, for instance, v+1 input samples. The effects resulting from the 
^ 10 transmission channel are modeled and represented by the h channel block 1012, where 
O the resultant, hk, denotes the original "long" channel impulse response. A: is a whole 

ftj number from 1 to Nq^ where Nq is the channel impulse response length of the channel h. 

The resultant signal hk is further distorted by additive disturbance v;^ (where ife is a real 
H number), which models the receiver noise and other signal interferences. This effect is 

.J 15 represented, for purposes of illustration, by the combination of hk and Vk at hypothetical 
adder 1002. The resultant is the combined, unfiltered signal yk that represents the signal 
as input at the receiver device, for instance a modem. The signal 3;^^ is received at the 
% TEQ filter 1014 for preprocessing as described further below. The signal yk is the 

H resultant vector of the current and previous received samples, N^y, and is passed through 

20 TEQ filter 1014, which operates to shorten the channel impulse response of yk and results 
in output signal vector Zk 1004. The TEQ vector \/ to be applied to the signal yk may be 
calculated using TEQ filter coefficients {wk}, and vice-versa, and is combined with the 
signal yk to produce the filtered output signal Zk 1004 having a shortened channel impulse 
response. 

25 Initially, the transmitted signal, Xk 1000, is assumed to be an uncorrelated and 

identically distributed wide sense stationary sequence having a zero mean and a squared 
variance of x^, Gx^. It should be understood that this assumption is not entirely accurate in 
that the existence of the cyclic prefix, which is the simple repetition of a portion of the 
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transmitted samples, causes some degree of correlation among the samples. However, 
generally, the separation between the repeated samples is large enough to neglect the 
correlation due to the prefix. 

Accordingly, the combined effects of the transmission filter, the receiver filter, 
and the distortion caused by the transmission channel are modeled by a linear time 
invariant system with the impulse response hk. As briefly described above, the receiver 
noise and other interferences are modeled by additive disturbance Vjt, which is assumed to 
be a wide sense stationary sequence having a zero mean and correlation coefficients 
represented by: 



where m is a real number and where is the noise autocorrelation sequence. 

In one manner, the TEQ filter vector is calculated using a number (designated 
as Nteq) of Finite Impulse Response (FIR) TEQ filter coefficients {wk}, where the 
coefficients are defined as: 



A goal of TEQ is to convert the channel impulse response hk into a target impulse 
response with effectively only a certain number of nonzero entries, designated as Nt^ 
Therefore, if the nonzero entries of the resulting channel are modeled with the sequence 
bk, where ^ is a whole number from 0 to Np-l , the target channel, can be formulated as: 



where d is the effective delay corresponding to the starting location of the non-zero 

segment of the impulse response. 



r^(m) = £(vjtv^+^) 



and 



w-[wi<r- wnteq f 
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For TEQ filtering in a DMT system, one optimal approach involves selecting 
{wk} such that the number of bits loaded per symbol, B, is maximized. This approach is 
equivalent to maximizing B as: 



where T is the set of transmit tones, 5/(w) is the desired signal energy at tone /, //(w) is the 
combined energy of ISI and ICI at tone i caused by the components of tn outside d<n< 
d+Nj-l, Ni(w) is the noise energy at tone /, and F is the effective gap which is a function 
of the constellation, coding gain, and the margin requirement. 

Unfortunately, the maximization of B with respect to w is a difficult, processor 
Q 10 intensive computation, and is impractical for most real world, adaptive implementations. 

nJ 

Therefore, some alternative near-optimal schemes have been developed. One such 
scheme includes maximizing the geometric Signal to Noise Ratio (SNRg) as: 

SNRg can be defined as the product of SNRs at each transmitted tone within the symbol. 
P 15 However, this process is also computationally complex and, consequently, is also 
unsuitable for practical applications. Alternative, comparatively lower complexity 
approaches such as that described in United States Patent No. 5,285,474, to Chow et al. 
have been developed. However, the reduction in computational complexity is achieved at 
an unacceptable cost of significant performance degradation or convergence instability in 
20 most real world applications. 

An alternative solution to the TEQ problem includes minimizing wall energy with 
a constant constraint on the power of the target impulse response which leads to a 
computationally demanding Singular Value Decomposition (SVD) based algorithm. 
Wall energy is defined as the residual impulse response energy outside of the region 
25 corresponding to the target impulse response. One manner explored for enabling such 
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wall energy mininiization includes utilizing a least square fit based approximation. 
Unfortunately, this solution has proven to result in unacceptable levels of performance 
degradation. 

An alternative solution to the above problem includes a conventional Minimum 
5 Mean Square Error (MMSE) approach, where the target impulse response (TIR) 
coefficients are assumed to be part of the desired final, filtered signal vector (the search 
vector). However, the squared norm constraint assumed in this approach to avoid an all- 
zeros solution leads to numerically difficult computations, such as an eigenvector 
computation to obtain TEQ coefficients. It is especially hard to implement direct 
10 adaptive training of the equalizer coefficients under this constraint. 

5:[ While theoretically it may be possible to achieve improved performance using a 

111 significantly longer TEQ filter, the computation needs, power and silicon area required to 
m implement such a long TEQ filter arrangement is not likely cost effective and largely 

impractical. Further, in practical applications, adaptive computation of TEQ filter 
^2 15 coefficients is a necessity. Still further, most digital signal processors (DSPs) in practical 

s . 

l^'. applications support only fixed point arithmetic. As a result, finding stable adaptive 

H algorithms with fixed point arithmetic to compute long TEQ filters may not be practical, 

p Therefore, there is a need in the art of digital ccwtnmunication systems for a practical 
method and system for effectively and efficiently shortening the channel impulse 
20 response hk to a filtered channel impulse response Zh 

The present invention overcomes the problems noted above, and provides 
additional advantages, by providing for a method of applying linear filtering to the 
received data set. The filtering is achieved by factoring the transmit data (xk) along a 
channel (hk) combined with additive disturbance/noise consideration (vk) to arrive at a 
25 combined unfiltered signal (yk\ and in consideration of a hypothetical sampling of the 
transmit data along a delay channel (d) to arrive at a target impulse response. From this, 
the present invention derives a Minimum Mean Square Linearly Constrained Time Delay 
Equalization (MLC-TEQ) mathematical model and applies the target impulse response to 
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achieve a shortened channel impulse response (zk)- The shortened channel impulse 
response defined by processing the unfiltered received data (yk) through the TEQ 
coefficients and TEQ vector is effectively compared with the hypothetical target 
channel impulse response (tn), which is derived from the delay channel as applied to (bk), 
to determine an error which is to be minimized by design. In this manner, the received 
data signal (yk) is converted into a TEQ filtered signal (zk) designed to approximate a 
desired target channel (tn). The invention shortens the overall channel response in a 
communication system such that the residual energy of the overall channel outside the 
target length is minimized. In addition, Inter Symbol Interference is minimized by the 
present invention. 

The block diagram of FIG. 10c illustrates the inventive MLC-TEQ filter approach 
that overcomes the shortcomings of earlier approaches and yields the aforementioned 
advantages. In particular, a transmit data signal, Xk 1000 is propagated along a channel 
1012 having distortive effects on the signal, as represented by the h channel block 1012, 
and has an effective channel impulse response represented as hk at 1020. To more fully 
represent the effects of various processes on the transmit signal, a disturbance vector Vk, 
which models the various receiver noises and interferences, at 1022 is applied to the 
chaimel impulse response to yield an unfiltered received data signal yk at 1024 that is 
received at an input of the TEQ filter 1014. By way of example, the received data signal 
yk could represent the input of a modem commonly used in communication systems. 
Unfiltered signal yk is the resultant vector of the current and previous Nw received 
samples, where Nw is defined as an integer representing the number of replicate received 
samples of transmitted data and represents TEQ length. 

The unfiltered yk signal is passed through Time Domain Equalization (TEQ) filter 
1014, which shortens the channel impulse response and generates a shortened channel 
impulse response signal Zk at 1026. As shown by the hatched lines, a hypothetical parallel 
branch of the transmitted data signal is shown as a delay function (d) 1032 applied to a 
detennined target vector bk to yield a target channel tn for modeling the shortened channel 
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impulse response to arrive at the TEQ filter coefficients. In effect, the TEQ filter 1014 
shortens the channel impulse response to minimize the mean square error ek, E{ek\ 
between the received signal and the hypothetical delay channel, thereby approximating 
the desired target channel. 

5 As shown in Fig. 10c, hypothetically, the delay channel 1032 generates a signal at 

1034 having a Target Impulse Response b 1036. The delay channel 1032 output is acted 
on by TIR b vector 1036 to produce the target impulse response channel tn is, in a sense, 
compared with the TEQ filter output vector Zk 1026 resulting in a data error signal eu at 
1004, which represents the ^h error signal (that is, the error signal that corresponds to the 

10 resultant data signal). In determining the desired TEQ filter coefficients, ek is minimized. 

In order to minimize the minimum mean least square error between the parallel 
communication channels (i.e., the impulse channel and the delay channel), it is necessary 
to calculate the TEQ filter coefficients used in generating the Equalizer Coefficient 
Vector w such that the resultant vector has minimum energy Inter Symbol Interference 
15 (ISI) and Inter Channel Interference (ICI). 

The MLC-TEQ approach of the present invention provides a family of solutions 
that are parameterized by variables d e {0,..., Nc + Njeq - Nt-2] and L e {0,..., Nj - 
1 }, It is further assumed that the channel impulse response {hi, i e {0,.. iVc - 1 } } and 
the noise autocorrelation sequence {rv(0, i ^ Z} are known or estimated. Each solution 
20 has a closed form expression in terms of d, {hi J e {0,..., A^c - 1 } }, {^v(0» ^ ^ Z}, TEQ 
length Nteq and the target impulse response length iVr. For fixed values of d and L, the 
MLC-TEQ filter determination steps are discussed in more detail below. 

The approach embraced by the TEQ filter determination process of the present 
invention poses the TEQ problem as one of minimizing the error differentiating the 
25 parallel branches of the h channel 1012 and hypothetical d channel 1032 of Fig. 10c to 
determine the desired TEQ coefficients {wk}, represented as: 

{Wk;kE{l Nteq }} 
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and Target Impulse Response coefficients {Z>k} as: 

{Z>k;ke{l,...^r}}. 

An additional constraint on w or b is used to avoid the trivial all zeros solution. In 
the present invention, that constraint is preferably on b, where one of the taps of b is 
constrained to be equal to a non-zero constant as = c, where L e { l...A/r}, ce% and c 
^ 0, C may be chosen as c = 1 without a loss of generality. This constraint allows the 
use of linear mean square minimization techniques, and the resulting algorithm is 
therefore more suitable for fast and efficient direct adaptive implementations. 

Due to the disturbance Vk and the finite length constraint on b and w, the error 

sequence {ek} has been found to generally not be equal to zero. Therefore, the error may 
be minimized in order to create effective filtering. In one manner, the mathematical 
method of the present invention, targets the minimization of the mean square of error 
In order to do so, the TEQ problem is preferably posed as to minimize Eiek') for w and b 
such that bi-l, where: 



'=[ 



b= 



and 



E(el) 

is the expected value (mean) of the square error, which is the minimization of mean 
square error by choosing vector variables w and b under the constraint that the Lth tap of 
vector b is constrained to be equal to 1. Minimization of E over w and the unconstrained 

components of ^ is a linear mean square minimization problem and the result for TEQ 
coefficient vector w is given as: 



-70- 



PATENT 

Attorney Docket No. 56162.000322 



Referring now to the flow chart of Fig. 11, one embodiment of a method for 
determining a set of TEQ coefficients to be used to generate a signal vector Zk in 
accordance with the present invention. Initially, in step 1100, the combination of the 
5 output of the h channel and the disturbance vector, jky is sampled. In step 1102, the 
hypothetical delay channel d is set. The impulse response length of the target channel (r„) 
is determined in step 1104. Optimal TEQ filter coefficients {wk} are determined in step 
1 106, such that the number of bits of data loaded per symbol {B) is maximized while the 
channel impulse response is minimized to provide an output with minimal noise and 
H 10 heightened data density. In step 1108, an equalizer coefficient vector is calculated 

ii : 

S based upon the TEQ filter coefficients defined in step 1106. In effect, the response Zk is 

Jjf compared, as in step 1110, with the target impulse response channel {tr^ vector b to 

rij generate optimal TEQ filter coefficients that result in minimizing the error of By 

% minimizing the error between the channels, the inventive method arrives at a set of 
f 15 coefficients Wk and vector to yield a shortened channel impulse response Zk that 

hi approximates the target channel impulse response. In a one embodiment, the Zk vector 

J!; preferably has a length Nt and is the target impulse response vector. 

2 Referring now to the flow chart of Fig. 12, one embodiment of the operation of 

step 1106 is described. In step 1202 the original long channel impulse response hk is 
20 determined. In step 1204, the noise autocorrelation sequence, Tv, is determined. In step 
1206, the delay value, rf, is fixed as the starting location of the non-zero segment of the 
target impulse response. Similarly, a value for L is chosen and fixed as the center of the 
range of possible values that L can take, where L is calculated as: 
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in step 1208, where L is an integer value used to set the target value of the length of the 
desired signal A^t-. With hk, r^^ d and L known or fixed at certain values, the error between 
the parallel branches Zk and bk may be minimized in step 1210. 

The specifics of step 1210 are discussed in more detail with reference to Fig. 13. 
In particular, Fig. 13 details an embodiment of a process used to calculate the TEQ filter 
{wk] and TIR filter {bk] coefficients. In order to minimize the error between Zk and bk, a 
MLC-TEQ mathematical method is constructed in step 1300. In order to avoid the trivial 
all-zeros solution, a finite length constraint is imposed onb ox w in step 1302. In the 
example of step 1302, one of the received signal channels, or "taps", of b is constrained 
as bi = c, where L e {0,...,A^r-l } and c e% and c 0. In this example, c is chosen to be 
c = 1; this is done without a loss of generality. 

A matrix F is constructed in step 1304 as: 



F= 



Ol,1 OlM-L-\ OL,NC+NTEQ-NT-d-l 

OnT-L-i ,1 InT'I^LNT-I^L Om-L^hNC+NTEQ-NT-dA 



A covariance matrix, Ry, is constructed in step 1306 as: 



^vNTEQ,NTEQ — 



rv(0) 
rv(l) 



rv(l) 
/V(0) 



MNteq-1) (Nteq-2) ... 



'•v(O) . 



A matrix H is constructed in step 1308 as: 



H = 



0 



0 Ao 



0 
0 
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Further, in step 1310, a matrix h is constructed for <i+L < Nteq as: 



h = [/? 



J + L 



h 



0 



0 



1 , NrEQ-d-L-l 



and for d+L > Nteq in step 1312 as: 




d + L 



h 



d + L - Nteq-1 



The Nteq FIR TEQ coefficients {wy^} are calculated in step 1314, as are the target 
impulse response coefficients {h} in step 1316. 

Returning now to Fig. 12, once the above values have been calculated, the various 

F, Rv, H and h and the coefficients {w^} are then used to calculate the Equalizer 
Coefficient vector in step 410 as: 



where: h is a function of the impulse response coefficient; Ry nteq,nteq is a noise 

2 

autocorrelation function; is the variance of H is a function of the channel 
impulse response; lis a function of the energy of the ISI and ICI interferences; F is an 
intermediate variable; and 



or the change from wq to w^teq -i . 

It should be understood that the choices of the parameters and L has a certain 
impact upon the performance of the TEQ filter with respect to a chosen application. 
Relating specifically to a DMT application, a reasonable capacity performance is 




J 
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obtained by choosing d as the location of the maximum value of the impulse response 

Further, it has been observed that for a typical scenario, there is a relatively wide 
region around dc where the performance of the filter does not change, where dc is the 
effective delay corresponding to the starting location of the non-zero segment of the 
impulse response. For the L parameter, a reasonable choice is given by the equation: 



L = 



where [xl is a function representing the smallest integer that is greater than or equal to 
Z, the number x. Therefore, L is chosen as the center of the range of possible values that L 

3 10 can take. 

>;l From the detailed description above, it will be clear to one skilled in the art that 

^ the present invention has many benefits. In particular, the fundamental advantages are 

^ two fold: it provides linear TEQ filtering with near-optimal performance, and it is 

structurally suitable for adaptive implementation on a real system. In particular, and in 
'2 15 contrast to the filtering methods currently used in the art, this method includes: the 
5 Target Impulse Response (TIR) vector is the combined response of the transmission 

channel and the TEQ filter is used as the search vector; the MLC-TEQ mathematical 
method constrains one element of the TIR vector to be equal to a constant and the 
element can be any element of the TIR vector; the error is computed using the difference 
20 between the TEQ output and the TIR filter output; the TEQ filter coefficients are 
computed by minimizing the Mean Square Error (MSE) criterion; and the ma±ematical 
method provides a family of TEQ vectors which are parametized by two quantities, the 
location of the constrained element of the TIR vector and the delay value used in the 
computation. 

25 In its Hadrian™ and Antoninus™ products, Virata Corporation of Santa Clara, 

California, extends the benefits of an integrated full-rate ADSL CPE chipset. The 
inventive concepts discussed above may be incorporated into an integrated ADSL CPE 
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processor and Analog Front End/Line Driver (AFE/LD) chipset, such as Virata 
Corporation's Hadrian™ and Antoninus™, which may be used in a wide variety of 
applications. According to one embodiment of the present invention, the inventive 
concepts related to echo cancellation functionality and TEQ functionality may be 
5 incorporated into the Hadrian and/or Antoninus products. 

FIGs. 14a and 14b are schematic diagrams of hardware architectures in which the 
inventive aspects of the present invention may be incorporated. FIG. 14a illustrates a 
block diagram of a chip architecture 1400 for Hadrian 1, which includes an integrated 
ADSL CPE processor 1402 and AFE/LD chipset 1440. The integrated ADSL CPE 

10 processor 1402 may be delivered as a multi-chip module (MCM) including Z3 chip 1412 
and Domitian Chip 1410. Z3 chip 1412 is full-rate ADSL Physical Layer (PHY) for CPE 
with on-chip DP (e.g., DSL PHY Processor) code storage memory. Domitian chip 1410 
is a high performance conmiunications processor for CPE with Ethernet PHY. Hadrian 
AFE/LD 1440 is an AFE with integrated line driver for full-rate ADSL CPE, as a part of 

15 the chip-set. The Hadrian AFE/LD 1440 is an integral part of the complete solution. The 
Hadrian Integrated ADSL CPE Processor 1402 is compatible with the Hadrian AFE/LD 
1440. 

The Hadrian- 1 chipset, as shown in FIG. 14a, is designed to be fully compliant 
with the full-rate ADSL ITU standard G.992.1 (G.dmt) and the splitterless ADSL ITU 

20 standard G.992.2 (G.lite). Hadrian- 1 fully supports Annex A, both Frequency Division 
Duplex (FDD) and Echo Cancelled (EC) modes, and Annex B of G.992.1 and Annex A 
and B of G.992.2. In addition, the initialization procedure is fiiUy compliant with 
G.994.1 (Rev. 2) (G.hs.bis). Other features may include being fully compliant with the 
splitterless ADSL ITU standard G.992.2 (G.lite) and supporting both Annex A and B. 

25 Hadrian 1 may also include a communications processor with Protocol Processor (PP) for 
protocol execution and Network Processor (NP) for data manipulation. Other features 
may include DSL PHY Processor (DP) for implementing modem firmware; Internal 
10/100 Ethernet PHY (disabled when Media Independent Interface (ME) interface is 
used); 10/100 Ethemet MAC with ME for connecting to either intemal 10/100 Ethernet 

-75- 



PATENT 

Attorney Docket No. 56162.000322 



saw:? 



PHY or external 10/100 multi-PHY; External Peripheral Bus (EPB); Synchronous 
Transfer Mode (STM): Serial frame data (e.g., High-Level Data Link Control (HDLC)); 
simultaneous STM and ATM support; bearer channels ASO, ASl, LSO, LSI for both 
STM and ATM; single and dual latency paths for both downstream and upstream; four 
5 framing modes, e.g., modes 0, 1,2, and 3; Network Timing Reference (NTR) transport 
support; dying gasp detection and handling; Universal Serial Bus (USB) 1.1 slave; 
integrated line driver and Voltage Controlled Crystal Oscillator (VCXO) for the AFE; 
general purpose input/output (GPIO); universal asynchronous receiver/transmitter 
(UART); In-Ckcuit Emulator (ICE); and/or Joint Test Action Group (JTAG). 

10 hi particular, standard compliance includes the power spectral density (PSD) of 

the transmit signals for all supported annex and corresponding modes; the use of 
:i frequency bins for transmit and receive signal including pilot bins; handshake procedure 
as specified in G.994.1; various initialization, training, and message signals; state 
II machine and timing sequence; DMT nK)dulation; framing modes; bearer channels and 
g; 15 dual latency paths; Embedded Overhead Control (EOC) channel; ADSL Overhead 
Control (AOC) Channel; forward error correction; inter-leaver operation; and ATM and 
IJ STM functionalities. Hadrian- 1 is interoperable with standard compUant central office 

g (CO) equipment, including various DSL Access Multiplexor (DSLAM) vendors. 

I FIG. 14b illustrates Hadrian 2, which includes an integrated ADSL CPE 

20 Processor 1452 and AFE/LD 1460 chipset, as shown by 1450. Hadrian 2 ADSL CPE 
processor 1452 may be delivered as a single-chip and cost-optimized device, which 
includes Z3 chip and Domitian chip. Z3 chip is a full-rate ADSL PHY for CPE without 
on-chip DP code storage memory. Domitian chip is a communications processor for CPE 
with Ethernet PHY and with SDRAM interface capable of handling various processors, 
25 such as Network Processor (NP), Protocol Processor (PP) and DSL PHY Processor (DP). 
Hadrian AFE/LD 1460 is an AFE wi± integrated Une driver for full-rate ADSL PHY for 
CPE and will also be included as part of a chip-set. Hadrian 2 may support the above 
listed features for Hadrian 1 in addition to STM: TDM H.110, ST-Bus. 
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The Domitian communications processor has been architected around the 
principle that there are two functions of a communications system which may use 
different architectures. These two functions are data manipulation and protocol 
execution. 

Data manipulation involves many word, byte and bit changes on small amounts of 
data. These operations have a deterministic character and are finished in a short time 
span. The code controlling these operations are small and compact, often written in 
assembly language, and involves a detailed knowledge of low level data protocols (Open 
Source Initiative (OSI) layer 1 to layer 3) and the available hardware. Data is processed 
in cells and packets of limited size. 

The protocol stack routines are multitasking, involving a large amount of memory 
and although time limits are set, these are not as strict as for the low level data 
manipulation routines. The data consists of frames which can be large (e.g., up to 64K). 

The Domitian architecture reflects these two types of function. The Protocol 
Processor (PP) has a 4 KB cache and is thus suited for running large programs. It has no 
direct access to the network data port but can exchange data with the Network Processor 
(NP) using shared Synchronous Dynamic Random-Access Memory (SDRAM). Large 
frames of data may be easily exchanged between the two processors using this method. 
The main PP tasks are to set up and close connections, get data from the NP, convert the 
data from one protocol type to another and pass them back to the NP again to be 
transferred to another network interface. For protocol support, it may also generate its 
own data frames for signalling, operating and maintenance. The PP may also 
communicate with, and control external hardware using the External Peripheral Bus 
(EPB). 

The Network Processor (NP) has a dedicated local program memory (e.g., 16 KB 
of SRAM) and may directly interact with network ports and associated hardware. A 
main tasks involves assembling incoming packets into frames for each of the network 
connections. It knows the status of each network port and responds to any incoming 
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packet/cell in time to prevent buffer overflow. While receiving packets/cells, it 
assembles them into frames, keeps track of errors, gathers statistical data, performs 
policing, extraction of Operation, Administration and Maintenance (0AM) cells, etc. 
When transmitting data, it segments frames into packets/cells, transmits at a specified rate 
5 by inserting idle cells, adds Cyclic Redundancy Check (CRC) bytes, inserts OAM cells, 
etc. 

FIG. 15 is a block diagram of a physical media dependent (PMD) layer 1500 of a 
ADSL CPE chip in which the inventive aspects of the present invention may be 
incorporated, FIG. 15 illustrates core building blocks for the PMD layer and basic signal 
10 flow. As shown in FIGs. 14a and 14b, PMD 1430 is implemented in ADSL CPE 
Processor 1402 and 1452. 

PMD layer may include a variety of functions, which may include a DSL PHY 
Processor (DP) for implementing modem firmware; a complete DMT engine including 
IFFT (e.g., 64-point, 128-point or multiples thereof) on the transmit path and a point Fast 
J: 15 Fourier Transform (FFT) (e.g., 512-point) on the receive path, Tx/Rx Bit Extract, Cyclic 
Prefix, and Constellation encoder/decoder; a Forward Error Correction (FEC) mechanism 
including the mandatory Reed Solomon (RS) encoder/decoder and optional Trellis 
m encoder/decoder; both fest and interleaved paths are supported along with the 

Interleaver/deinterleaver; Training Signals Generator and Sync Symbol Generator; 
20 Course/Fine Gain Scaling and Tx/Rx Clip Scaling; Programmable Transmit Filter and 
Receive Noise Reduction Filter; Digital Automatic Gain Control (AGC); Adaptive Time 
domain Equalization (TEQ); Adaptive Frequency domain Equalization (FEQ); Tone 
ordering and Bit Swap; Clock Recovery Loop (CRL); Sigma-Delta Interpolator and 
Decimator for Digital to Analog Converter (DAC) and Analog to Digital Converter 
25 (ADC); and Adaptive Echo Cancellation Filter. 

As shown in FIG. 15, adaptive echo canceller 1510 may embody the inventive 
aspects related to echo cancellation of the present invention discussed above. In addition. 
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adaptive TEQ 1520 may embody the inventive aspects of the TEQ of the present 
invention discussed above. 

FIG. 16 is a block diagram of a transmission convergence (TC) layer 1600 of a 
ADSL CPE chip in which the inventive aspects of the present invention may be 
5 incorporated. As shown in FIGs. 14a and 14b, TC 1420 is implemented in ADSL CPE 
Processor 1402 and 1452. The TC layer may support a variety of features, such as ATM- 
TC; STM-TC; dual-bearer; simultaneous ATM and STM transport; Bearer Channel 
Switch; UTOPIA 2 interface; ADSL Overhead Control/Embedded Overhead Control 
(AOC/EOC) Buffers; Tx/Rx CRC; and Scrambler/Descrambler. 

10 FIG. 17 is a block diagram of an analog front end device in which the inventive 

aspects of the present invention may be incorporated. Hadrian AFE/LD, as shown by 
1440 and 1460 in FIGs. 14a and 14b, with integrated line driver is an integral part of the 
complete solution. To achieve enhanced performance the Hadrian AFE/LD and its 

Hadrian Integrated ADSL CPE Processor counterpart are designed with tightly coupled 
15 architecture. AFE features may include an integrated line driver; an integrated VCXO; 
16-bit DAC with 14.5 effective bits; 16-bit ADC with 13.5 effective bits; industrial 
temperature range: -40C ~ 85C; and external filters. 

FIGs. 18a- 18c are a block diagram of applications in which the inventive aspects 
of the present invention may be incorporated. FIG. 18a illustrates an Ethernet 
20 Router/Bridge 1800 which incorporates the Hadrian product. FIG. 18b illustrates a USB 
Modem and USB- Attached Gateway 1820 which incorporates the Hadrian product. FIG. 
18c illustrates an Integrated Access Device (IAD) 1840 which incorporates the Hadrian 
product. 

Virata's Antoninus is a full-rate ADSL CPE PHY device comprised of two chips, 
25 which includes an ADSL CPE Processor called Z3 (as discussed above in connection 
with Hadrian 1) and a AFE/LD chipset, as illustrated in FIG. 17. The Z3 chip is a full- 
rate ADSL PHY for CPE with on-chip DP (e.g., DSL PHY processor) code storage 
memory. The ADSL CPE processor Z3 is also used in Hadrian 1, discussed above. 
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Similarly, the inventive aspects related to the echo canceller and the TEQ may be 
implemented in the Antoninun product. Antoninus AFE/LD is an AFE with integrated 
line driver for full-rate ADSL CPE, which is also included as part of the chip-set. 

While the foregoing description includes many details and specificities, it is to be 
understood that these have been included for purposes of explanation only, and are not to 
be interpreted as limitations of the present invention. Many modifications to the 
embodiments described above can be made without departing from the spirit and scope of 
the invention. 

The present invention is not to be limited in scope by the specific embodiments 
described herein, hideed, various modifications of the present invention, in addition to 
those described herein, will be apparent to those of ordinary skill in the art from the 
foregoing description and accompanying drawings. Thus, such modifications are 
intended to fall within the scope of the following appended claims. Further, although the 
present invention has been described herein in the context of a particular implementation 
in a particular environment for a particular purpose, those of ordinary skill in the art will 
recognize that its usefulness is not limited thereto and that the present invention can be 
beneficially implemented in any number of environments for any number of purposes. 
Accordingly, the claims set forth below should be construed in view of the full breath and 
spirit of ±e present invention as disclosed herein. 
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